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ABSTRACT 


This  survey  reveiws  the  analytical  literature  pertinent  to  contained 
rotating  flows.  A  related  survey  of  the  experimental  literature  Is  re¬ 
ported  In  TM  67-2.  The  main  emphasis  In  this  report  Is  placed  upon  posslbl 
applications  to  the  study  of  the  fluid  mechanics  In  rotating  rocket  motors. 
Since  this  area  is  only  now  receiving  analytical  attention  It  was  necessary 
to  draw  heavily  upon  the  classical  analytical  studies  of  vortex  tubes  and 
swirling  flows  in  cylinders  and  converging  nozzles.  A  great  deal  of  the 
classical  fluid  mechanics  in  this  area  is  incompressible  In  nature,  there¬ 
fore,  must  be  considered  qualitative  In  relation  to  actual  rocket  motors. 
Conditions  existing  In  the  nozzle  and  chamber  are  reviewed.  In  addition 
other  topics  such  as,  end  wall  interactions  and  boundary  layers  are  re¬ 
ported  upon.  The  analytical  models  for  the  determination  of  the  effects 
of  acceleration  on  the  burning  rate  of  solid  propellants  which  have  been 
reported  in  two  recent  papers  are  reveiwed. 


1.0  INTRODUCTION 


The  material  presented  In  this  literature  survey  will  consist  of 
analytical  attempts  to  describe  the  difficult  problem  of  contained 
rotating  flows.  The  main  emphasis  is  placed  upon  the  theory  of  vortex 
flows  In  rocket  motors.  As  yet,  there  Is  not  a  great  deal  of  analytical 
work  specifically  oriented  towards  vortex  flows  in  rocket  motors.  How¬ 
ever,  there  has  been  a  great  deal  of  interest  in  vortex  flows  In  general 
over  the  years.  Therefore,  many  of  the  articles  reviewed  must  be  con¬ 
sidered  supplementary  In  nature,  yet  most  of  the  fundamental  points 
pertinent  to  rocket  vortices  are  covered. 

The  discovery  of  the  Ranque-Hllsch  vortex  tube  stimulated  a  renewal 
of  Interest  In  vortex  flow  theory.  Many  of  the  theories  developed  to 
understand  this  device  have  applicability  to  the  conditions  existing  In 
the  chamber  of  rotating  solid  propellant  rocket  motors. 

The  study  of  rotating  disks  provides  an  Insight  Into  the  problem  of 
the  Interaction  of  the  end  walls  with  chamber  flow  field.  The  long  famlll 
bathtub  vortex  In  which  a  vortex  is  generated  and  flows  toward  a  sink 
prompted  early  interest  In  vortex  flow  in  nozzles.  Much  of  the  early  work 
was  for  an  Incompressible  perfect  fluid.  It  was  noted  that  the  flow  rate 
In  such  a  nozzle  was  less  than  that  for  a  non-rotating  nozzle.  Later, 
the  addition  of  compressibility  effects  made  It  necessary  to  consider  a 
void  region  near  the  axis,  since  the  potential  vortex  possess  the  tendency 
toward  Infinite  centerline  velocities.  The  addition  of  viscous  effects 
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further  complicated  the  theory  In  that  a  solid  body  type  vortex  Is 
necessarily  established  at  the  center  line  so  that  the  free  boundary 
conditions  at  the  axis  are  satisfied. 

Recently,  the  possibility  of  the  existence  of  regions  of  reversed 
axial  flow  has  been  Introduced  to  explain  some  apparently  anomolous  ex¬ 
perimental  results.  Flow  reversal  has  been  attributed  to  both  viscous 
effects  such  as  end  wall  boundary  layer  Interactions  with  the  primary 
flow  field,  and,  dynamic  effects  wherein  finite  transitions  occur  between 
two  conjugate  flows.  The  latter  of  these  Is  referred  to  as  vortex  break¬ 
down  and  Is  usually  associated  with  a  stagnation  region  on  the  axis  with 
locally  reversed  flow. 

The  literature  survey  also  Includes  two  published  theories  of  burning 
rate  augmentation  In  spinning  rocket  motors.  The  fact  that  the  propellant 
Is  burning  differently  from  the  non-rotating  case  Is  important  In  establish¬ 
ing  the  Initial  conditions  for  any  detailed  study  of  the  flow. 

Finally,  two  articles  are  presented  In  which  the  entire  rocket  motor, 
the  chamber,  the  end  walls,  and  the  nozzle  are  considered  to  determine 
the  ability  of  the  vortex  to  effect  the  overall  performance  of  the  motor. 


2.0  COMPRESSIBLE  VORTEX  FLOWS  IN  NOZZLES 


2.1  Introduction 

The  flow  of  a  vortex  In  a  nozzle  has  been  studied  for  a  great  many 
years.  Historically  the  first  theories  were  presented  for  the  Invlscld, 
trrotational  flow  of  liquids.  It  was  apparent  from  these  calculations 
that  the  mass  flow  rate  in  a  nozzle  In  which  the  fluid  was  rotating  was 
not  as  great  as  the  flow  rate  In  a  non-rotating  case  with  similar  chamber 
conditions.  The  formation  of  an  air  core  along  the  axis  of  the  liquid 
also  served  to  reduce  the  mass  flow  rate.  Later  when  compressible  flow 
theory  for  gases  was  used  an  analogous  situation  developed  In  which  the 
mass  flux  at  the  throat  was  reduced  due  to  the  rotating  component. 

Slmlllar  to  the  Incompressible  case  a  void  region  near  the  axis  was  pre¬ 
dicted  since  the  Irrotational  vortex  predicts  infinite  tangential  velocities 
at  the  center  line.  However,  conservation  of  energy  dictates  that  when 
the  tangential  velocity  reaches  the  maximum  isentropic  speed  further  ex¬ 
pansion  Is  Impossible  since  the  pressure  is  zero. 

Recently,  in  connection  with  end-burning  rocket  motors,  the  flow 
of  a  rotational  vortex  has  been  studied.  These  studies  indicate  a  slmiliar 
decrease  in  the  expected  mass  flow  rate  for  given  chamber  conditions. 
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2.2  Irrotational  Vortices 

2.2.1  Theory  of  Binnie  and  Hooker 

In  a  paper  In  1937,  Binnie  and  Hooker  (1)  discuss  the  radial  and 
spiral  flow  of  a  compressible  fluid  in  which  a  velocity  potential  is 
assumed  to  be  definable, 

$  =  A0  +  B  f(r)  (2.1) 

Since  a  velocity  potential  Is  assumed,  the  flow  field  must  be 
considered  to  be  irrotational  and  homentropic,  Therefore,  the  proper 
expressions  for  the  tangential  velocity  and  radial  velocity  are  given  by, 

v  ■  *  If  •  u  -  & 

It  was  shown  that  for  radial  flow  of  a  swirling  gas,  the  radial 
velocity  must  be  less  than  In  the  case  of  no  swirl.  In  addition,  the 
radial  velocity  does  not  reach  the  sonic  condition  until  the  minimim 
cross  sectional  area  Is  reached  where  the  combined  velocity  of  the  gas 
is  greater  than  the  sonic  velocity.  The  minimum  area  corresponds  to  the 
point  where  the  area  must  begin  Increasing  if  the  velocities  are  to  con¬ 
tinue  increasing. 


2.2.2  Theory  of  Binnie 

In  1949,  Binnie  (2)  again  reported  on  the  invlscid,  adiabatic  flow 
In  an  axlsymmetric  nozzle.  In  addition  to  other  choking  phenomena,  the 
flow  of  a  spiraling  gas  In  a  nozzle  was  discussed.  As  in  previous  work, 


5 


■3f 


4 


■:* 


4 


the  tangential  velocity  was  assumed  to  be  given  by, 

V  =  c/r  (2.3) 

so  that  the  flow  was  irrotational .  It  was  also  assigned  that  the  radial 
velocity  was  everywhere  small  so  that  it  could  be  neglected. 

Employing  the  Euler  equations  for  axisymmetric,  inviscid  flow  and 
neglecting  radial  velocities,  Binnie  derived  the  equation  for  mass  flow 
in  the  nozzle. 

i  =  fw  W0o  [1  -  41  t**2  *  2,rdr  (2.4) 

ro  ao 

where  rQ  =  the  radius  of  the  inner  void  region 
r  =  the  radius  of  the  nozzle  wall 

W 

aQ  =  the  stagnation  accoustic  speed 

For  analytical  convenience,  Binnie  chose  y  =  2.0  to  study  the  character 
of  the  equations  at  the  throat  condition.  Although  not  stated,  it  was 
assumed  that  the  axial  velocity  was  a  constant  at  every  cross  section  in 
the  nozzle.  Therefore,  he  was  able  to  integrate  the  expression  for  mass 
flow  rate.  Then  by  differentiating  that  result  with  respect  to  the  axial 
dimension  and  setting  the  result  equal  to  zero,  the  conditions  at  the 
"throat"  were  found.  It  was  shown  that  when  c  -  o,  the  equation  for  the 
axial  velocity  reduces  to  the  local  speed  of  sound  as  in  the  case  for  one¬ 
dimensional  flow.  Binnie  further  illustrated  that  when  c  f  o,  the  axial 
velocity  must  be  some  value  less  than  the  local  speed  of  sound.  This 
value  of  axial  velocity  was  related  to  the  value  of  the  parameter,  which 
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Is  the  angular  momentian  of  a  streamline. 

Bfnnle  presents  figures  for  the  reduced  axial  velocity  In  the  throat 
as  a  function  of  the  ratio  of  the  nozzle  to  void  region  radii  and  the 
angular  momentum  of  a  streamline.  The  figures  were  presented  for  y  =  2  and 
rs  1.4. 

In  conclusion,  this  analysis  was  presented  for  an  inviscid,  Irrotational , 
adiabatic  flow  in  which  radial  velocities  were  neglected  and  axial  velocities 
were  assuned  uniform  over  each  cross  section. 

2.2.3  Theory  of  Mager 

In  1961,  Mager  (3)  made  further  contributions  to  the  theory  of  swirling 
flow  In  a  nozzle.  As  in  the  earlier  work  by  Blnnie,  Mager  assumed  an  Ir¬ 
rotational,  Inviscid,  adiabatic  flow  in  which  radial  velocities  are  neglected 
and  the  axial  velocity  is  assumed  uniform  In  the  radial  direction.  A 
schematic  of  the  physical  model  Is  presented  In  Figure  (1). 

In  his  work,  Mager  assumed  a  velocity  potential  function  of  the  follow¬ 
ing  form: 

*  *  »(z)  w(z)’9  +  I  w(z)dz  <2-4> 

so  that 

u  =  Dr  -  0;  W  =  (RW)ze  +  w;  v  =  ^  60  ,  (2.5) 

and  the  total  velocity  is  then  given  by 

q2  3  (c2/r2  +  w2)  . 


(2.6) 
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The  mass  flow  rate  through  the  nozzle  is  specified  In  the  same  manner 
as  that  determined  by  Binnie 

m  =  fr  wpQ  [1  -  ^  (c2/r2  +  w2)]1/y_1  2irrdr  =  K2  (2.7) 

iro  2a 

since, 


V  *  c/r  . 


(2.3) 


It  can  be  seen  that  as  r  approaches  zero,  the  tangential  velocity  tends 
to  Infinity.  However,  the  velocity  of  the  gas  is  physically  limited  to 
the  maximum  Isentropic  velocity  corresponding  to  zero  temperature.  Thus, 
an  expression  can  be  derived  for  the  radius  of  the  void  region 


r 


o 


2 


(2.8) 


The  equation  for  mass  flow  can  then  be  Integrated  for  the  selected  values 
of  y  which  can  be  expressed  In  the  following  form. 


Y 


(2.9) 


Expressions  for  mass  flow  rate  and  specific  impulse  were  then  derived 
as  a  function  of  the  parameter 


a  s  (c/R  )  ^  Y-l  • 
ao  2 


(2.10) 


Charts  were  presented  illustrating  the  effect  of  the  value  of  a 
at  the  "throat"  upon  the  mass  flow  and  thrust  as  compared  to  one-dimensional 


3  |3 


9 


flow.  It  was  shown  that  the  greater  the  contraction  ratio  of  the  nozzle, 
the  greater  the  reduction  of  the  mass  flow  and  specific  Impulse  (Figure  1). 

2.3  Rotational  Vortices 

2.3.1  Theory  of  Crocco  and  Vaszonyl 

It  was  shown  by  Crocco  and  Vaszonyl  (4,5,6),  that  when  gradients  exist 
In  either  the  stagnation  temperature  or  the  stagnation  pressure,  the  flow 
of  the  fluid  is  in  general  rotational.  The  generalized  unsteady  Crocco- 
Vaszonyi  equation  states 

TvS  +  q  x  V  x  q  =  VH  + 

O  dt 

and,  for  steady  state, 

q  x  V  x  q  =  q  x  iL  =  VHq  -  TVS 
There  are  three  conditions  under  which  q  x  to  =  0  .  First, 
q  =  0  which  Is  a  trivial  case. 

Second, 

q  x  q  =  u  =  0  which  is  the  case  of  an  Irrotatlonal  flow. 

Third, 

q  ||  to  -  which  is  a  special  case  of  rotational  vortex  motion  in 
which  the  vortex  lines  are  everywhere  parallel  to  the  streamlines.  This 
type  of  motion  Is  called  Beltrami  Helical  motion  and  can  be  derived  from 
a  uniformly  translated,  two-dimensional  vortex  In  a  constant  area  section. 


(2.11) 


(2.12) 


(2.13) 
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Thus,  fn  general,  for  the  flow  of  a  gas  which  has  gradients  In  the 
stagnation  properties, 


vs  =  cp  v  In 


^  *1 


0  Y 


(2.14) 


Thus,  the  entropy  Is  not  constant  in  the  direction  normal  to  the  stream¬ 
lines,  and  the  flow  Is  rotational.  Substituting  eqn.  2.14  into  eqn.  2.12 
gives  the  following  result, 


qxuj=cpvT0-Tv1n  ( - 2__)  f  o 

Po  *7“ 


(2.15) 


except  when,  TQ  and  P0  are  constants.  Therefore  If  gradients  exist  In 

either  T  ,  or  P„  the  flow  field  must  be  considered  rotational, 
o  o 

A  further  consequence  of  variations  of  stagnation  properties  is  that 
the  well  known  law  of  Kelvin  In  which  the  circulation  for  a  barotropic 
fluid  Is  constant  for  all  time,  Is  no  longer  applicable.  A  barotropic 
flow  exists  when  there  is  a  unique  relationship  between  the  pressure  and 
density.  The  two  most  Important  examples  are  Incompressible  flow  and 
isentroplc  flow.  In  a  flow  which  has  variations  In  the  stagnation  pro¬ 
perties,  there  no  longer  exists  a  unique  relationship  between  the  pressure 
and  density  normal  to  the  streamlines.  However,  along  a  streamline,  the 
flow  may  still  be  isentroplc. 

Thus,  the  circulation  r,  which  Is  given  by  the  expression 


r  =  o  V  *  dl 


(2.16) 


is  no  longer  a  constant  in  each  plane. 


2.3.2  Theory  of  Bastress 


The  first  study  of  the  flow  of  a  rotational  vortex  In  a  nozzle  was 
presented  by  Bastress  (7)  In  1965.  This  analysis  was  followed  by  several 
others  which  are  essentially  similar.  However,  each  exhibits  certain  re¬ 
finements.  Since  the  stagnation  properties  are  In  general  not  constant 
normal  to  the  streamlines,  numerical  calculations  associated  with  these 
analyses  rely  upon  average  properties  within  stream  shells.  A  stream 
shell  Is  defined  as  an  annulus  through  which  a  given  amount  of  mass, 
energy  and  momentum  are  flowing  and  across  which  the  stagnation  properties 
are  averaged,  as  Illustrated  In  Fig.  2. 

The  following  set  of  assumptions  Is  common  to  all  the  analyses  pre¬ 
sented  In  this  section. 

1.  Steady  flow 

2.  Invlscid  flow 

3.  Adiabatic  flow 

4.  Perfect  gases 

5.  Negligible  radial  velocities 

It  assumed  either  directly  or  by  implication  that  a  solid  body  vortex 
persists  at  every  cross  section  of  the  nozzle.  In  addition,  all  analyses 
are  basically  one-dimensional  In  that  only  radial  variations  of  the  pro¬ 
perties  of  the  flow  are  allowed.  However,  both  tangential  and  axial 
velocities  are  considered.  A  schematic  of  the  physical  model  Is  given  In 
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In  the  work  by  Bastress,  the  following  equations  are  considered. 


State 

P  =  pRt 

(2.17) 

Energy 

2  2 

c  T  =  c  t  +  w  t  v 
pop  2 

(2.18) 

Angular  momentum 

“oro2  =  V  *  r 

(2.19) 

Entropy 

VS  =  0 

(2.20) 

In  addition  to  the  general  assumptions  listed  above,  Bastress  also 
assumes, 

1.  There  exists  uniform  contraction  ratios  for  the  stream  lines. 

2.  The  flow  Is  hometropic. 

3.  The  radial  pressure  gradient  Is  zero. 

The  first  assumption  Implies  solid  body  rotation  at  each  plane.  This 
can  be  shown  by  considering  the  conservation  of  angular  momentum  for  two 
streamlines  within  the  flow. 

On  the  streamline,  conservation  of  angular  momentum  gives. 


r  2  -  r  2 
oj,r  .  -  (u-r. 

o  OJ  J  J 


(2.21) 


On  the  j  +  1th  streamline,  conservation  of  angular  momentum  gives, 


2  2 
°oroj+l  '  U)j+lrj+l  * 


(2.22) 


Thus , 


oj.  =  uj  0  JO 

J  o  — s—  » 


“j  +  1 


Oj  +  1 

0  r  2 
rj+l 


(2.23) 
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Since  the  contraction  ratios  of  each  streamline  are  assumed  constant, 
_  2 


“o  "  “j 


r.2 

=  u.+i  ^  •  =  a  constant  with  radius. 
roj  roj+l 


(2.24) 


The  second  assumption  leads  to  the  conclusion  that  the  stagnation 
temperatures  and  pressures  should  be  considered  constant  normal  to  the 
streamlines.  However,  this  violates  the  Crocco-Vaszonyl  equation  dis¬ 
cussed  earlier.  This  contradiction  arises  because  a  solid  body  vortex 
is  In  general  rotational  even  when  the  gradient  of  the  vortlclty  vector 
Is  zero,  as  Is  found  In  a  Beltrami  Helical  motion. 

The  third  assumption  leads  to  the  conclusion  that  centrifugal  forces 
are  negligible  and  that  the  thermodynamic  variables  may  not  vary  with 
radius.  An  additional  consequence  of  the  second  and  third  assumption  Is 
that  the  total  velocity  must  be  a  constant  In  each  plane  normal  to  the 
axis. 

Based  upon  these  assumptions  and  the  listed  equations,  Bastress 
developed  an  equation  for  the  change  In  effective  throat  area  with 
angular  momentum  of  the  flow.  The  axial  velocity  at  the  throat  Is  given 
by, 


w*  =  (a*2  -  $*2r2)1/2 


(2.25) 


where, 

<f>  =  (v  -  u)Qr) 

a*  =  the  accoustlc  speed  at  the  throat 
o>0  =  the  Initial  spin  rate 
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■M 


M 


Ml 


•* 


ta 


A  * 
Aeff 

~W~ 


1 


rD*/2  .*2  2  ,  ,/9 

2irr(l  -  ^-£-)dr  =  |  2  [1  -  (l-e)3/2] 
J  o  a* 


where. 


e  *  **  D*/2a*  =  (  QD*/2a*)  [(DQ/D*)2  -  1]  . 


(2.26) 


The  properties  of  Eqn.  (2.26)  are  such  that  as  the  spin  rate,  a)Q,  goes 
to  zero  the  ratio  of  the  areas  approach  unity.  The  parameter  e  has  a 
maximun  value  of  unity.  At  this  point,  the  effective  area  reduction  Is 
two-thirds.  When  e  approaches  unity,  the  axial  velocity  at  the  wall  approaches 
zero.  Thus,  there  Is  a  limitation  to  the  analysis  In  that  the  initial  spin 
rate  cannot  exceed  a  certain  value  because  of  the  Initial  assumption  of  uniform 
contraction  ratios. 

Bastress  then  relates  the  change  In  effective  throat  area  to  a  change 
In  chamber  pressure  using  conventional  Internal  ballistic  relations.  Thus, 

It  Is  concluded  that  rotation  can  reduce  the  effective  throat  area,  and 
thereby  Increase  the  chamber  pressure  and  burning  rate  of  the  propellant  as 
illustrated  in  Fig.  3. 


2.3.4  Theory  of  Manda 

An  additional  study  related  to  the  flow  In  a  nozzle  generated  by  an 
end-burning  rocket  motor  Is  presented  by  Manda  (8).  In  this  work,  certain 
Improvements  are  made  upon  the  theory  of  Bastress  by  Including  the  radial 
momentum  equation,  thereby  eliminating  the  necessity  for  the  assumption  of 
constant  total  velocity  and  the  Implied  assumption  of  uniform  thermodynamic 
properties  across  a  plane.  The  assumption  of  uniform  contraction  ratio. 
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however.  Is  retained. 

The  pertinent  equations  of  this  study  are  given  below. 


State 

Energy 

Angular  Momentum 
Radial  momentum 
Entropy 


P  »  PRt 


(2.27) 


2  2 

cpT0  s  c  t  +  —  2  v-  =  const,  on  streamlines  (2.28) 


Vr  =  r  -  const,  on  streamlines 

1  aP  = 
p  ar  r 


(2.28a) 

(2.29) 


VS  =  0 


It  might  be  pointed  out  that  the  lack  of  a  continuity  equation  is 
replaced  by  the  assumption  of  uniform  contraction  ratios  which  implies 
that  all  the  flow  in  the  chamber  is  conserved  in  the  stream  shells.  There¬ 
fore,  the  location  of  the  streamlines  may  be  determined  from  the  contraction 
of  the  nozzle  wall. 

By  combining  the  equations  for  energy  and  angular  momentun  In  a 
stream  shell,  along  with  an  isentropic  expansion  from  the  chamber  conditions 
(T0,  PQ),  Manda  obtains  the  following  expression. 

2gyRT„  T  o  o  m  r  n 

— TTf-  (f  -  1)  +  W2  -  WQ2  +  (r0u,)2  [02  -  1]  =  0  (2.31) 

'  o 


Since  the  axial  velocities  emlnatlng  from  the  propellant  surface  are 

2 

small,  the  axial  velocity  WQ  is  subsequently  neglected.  Therefore,  an 
expression  for  the  axial  velocity  can  be  obtained. 


2gRyT 

~~y-T 


Rn  ,  R„  ,  , 

(rt)  Up2)  -  1]  r2 


(2.32) 
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w  here. 


R0/R  =  rQ/r 


Manda  next  employs  the  radial  momentum  equation  which  Is  Integrated 
assuming  Isentroplc  relations  between  the  thermodynamic  variables. 


1  BP  _  r 
p  Br  "  r 


(2.29) 


Although  not  stated,  this  Integration  then  implies  the  existence 
of  a  homentroplc  flow  field.  However,  as  In  Bastress1  analysis,  this  Is 
a  contradiction  to  the  fact  that  a  solid  body  vortex  Is  rotational. 

Thus,  the  following  expression  Is  obtained  at  the  nozzle  throat, 

p 

w*  =  a*2  -  u,02  (Rq/R*)  [2  (^)2  -  1]  r2  (2.33) 


where  R*  =  radius  of  the  geometrical  throat 


R0  =  chamber  radius 


=  the  Initial  motor  spin  rate. 

Substituting  this  expression  Into  the  equation  for  the  axial  velocity 
evaluated  at  the  throat  gives. 


Ill 

(P/P0)  Y  = 


(I  )  =  _£_  +  Id 
V  Y+l  2grRT0 


(2.34) 


Manda  then  applies  the  result  to  an  equation  for  the  mass  flow  rate 
at  the  throat  so  that  the  effect  of  the  initial  angular  momentum  upon 
the  mass  flow  rate  can  be  determined. 
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«*1 2 3 4 5  -  a.2  -  »02  {jjj)2  [2  t^-)2  -  1]  r2  (2.35) 

In  conclusion,  Manda  finds  that  the  mass  flow  rate  can  be  significantly 
reduced  for  large  values  of  Initial  angular  momentum.  However,  the  effect 
Is  reduced  for  a  high  flame  temperatures  of  the  propellant  gases.  Manda' s 
results  are  presented  in  Fig.  4. 

Summary 

Recently,  King  (9)  has  published  a  discussion  of  the  two  analyses  of 
end  burning  motors  presented  by  Bastress  and  Manda.  In  this  review,  some 
of  the  contradictory  statements  and  results  of  these  analyses  are  discussed. 
King  concludes  that  the  basic  error  in  the  two  analyses  is  the  assumption 
of  uniform  contraction  ratio.  This  assumption  leads  to  an  overspecification 
of  the  problem.  This  can  easily  be  seen  by  summarizing  the  equations 
employed  by  each  author. 


Bastress 


1.  State 

2.  Energy 

3.  =  0 

a  r 

4.  Angular  Momentum 

5.  Entropy,  vs  -  o 


Unknowns 

p,  p,  t,  u,  v  =  f ( r ) 


FIGURt  3  THE  RESULTS  OF  THE  THEORY  OF  BASTRESS(4). 
(Taken  from  reference  (4)  ,  Bastress) 


PROPERTY  DISTRIBUTION  AT  THE  THROAT 


BURNING  RATE  EXPONENT  INFLUENCE 
ON  COMBUSTION  PRESSURE 


FIGURE  4  THE  RESULTS  OF  THE  THEORY  OF  MANDA(5). 
(Taken  from  reference  (5)  ,  Manda) 
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Manda 

Equations  Unknowns 

1.  State  p,  p,  t»  u,  v  =  f(r) 

2.  Energy 

3.  Radial  Momentum 

4.  Angular  Momentim 

5.  Entropy,  vs  =  o 

In  addition  to  the  equations  listed,  each  analysis  employs  the 
assumption  of  uniform  contraction  ratios.  As  stated  earlier,  this 
assumption  Is  equivalent  to  solid  body  rotation  at  every  plane.  Mathe¬ 
matically,  this  means  that  the  dependent  unknowns  are  reduced  to  four. 
Thus,  it  can  be  seen  that  there  exists  one  too  many  equations  for  the 
number  of  dependent  unknowns.  Thus,  each  analysis  is  over-specified. 

Further  evidence  of  over  specification  is  that  each  analysis  employs 
a  different  set  of  governing  equations.  A  consistent  set  of  equations 
are  given  below, 

1 .  State 

2.  Energy 

3,  Angular  Momentum 

4,  Radial  Momentum 

5:  Mass  flow 

6.  Entropy 

Thus,  it  can  be  seen  that  there  are  six  equations  and  six  unknowns 
which  are  a  function  of  radius.  Therefore,  no  assumptions  are  necessary 
regarding  the  nature  of  the  tangential  velocity  (68). 


In  addition  to  the  problem  of  overspecification,  there  is  a  contra¬ 
diction  in  the  treatment  of  the  entropy  of  the  flow.  Each  of  the  authors 
assumes  a  homentroplc  flow  field  either  by  implication  or  direct  statement. 
However,  It  was  shown  in  an  earlier  section  that  a  solid  body  vortex  Is 
rotational  In  nature.  Hence,  a  variation  in  the  entropy  normal  to  the 
streamlines  is  necessary.  An  additional  consequence  of  the  rotational 
nature  of  the  flow  Is  that  the  stagnation  properties  must  vary  normal  to 
the  streamlines. 

Thus  a  complete  approach  to  the  problem  must  Include  entropy  gradients 
and  the  exlstance  of  stagnation  property  variations  normal  to  the  stream- 
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3.  VORTEX  FLOWS  IN  ROCKET  MOTOR  CHAMBERS 
3.1  Introduction 

The  properties  of  the  flow  In  a  rocket  motor  chamber  are  very  Important 
to  the  overall  study  of  vortices  in  rotating  rocket  motors,  since  these 
properties  serve  as  initial  conditions  for  the  flow  in  the  nozzle.  In  addition, 
the  flow  field  in  the  chamber  provides  information  as  to  the  effects  of  the 
vortex  on  propellant  burning. 

In  general  the  vortex  flow  in  an  end-burning  motor  is  easily  described, 
since  in  that  case  the  flow  is  generated  with  a  solid  body  vortex.  Under 
ideal  conditions  this  solid  body  vortex  would  travel  down  the  chamber  and 
enter  the  nozzle  without  deformation,,  The  action  of  viscosity  would  not 
necessarily  disrupt  the  solid  body  vortex  since,  when  the  walls  are  rotating, 
there  is  no  tangential  velocity  boundary  layer  at  the  walls.  Although  the 
solid  body  vortex  Is  considered  to  be  an  equilibrium  flow.  It  is  not  necessarily 
stable  to  disturbances  (53).  Recently,  Hall  (10)  presented  indications  that 
Instabilities  may  develop  If  gradients  in  the  axial  velocity  are  present. 

In  the  case  of  an  internal-burning  cylindrical ly  perforated  core,  the 
theory  is  considerably  more  complicated.  In  the  ideal,  inviscid  case,  a  free 
vortex  pattern  would  be  established  since  every  particle  of  gass  is  generated 
at  the  same  radial  location  and  thus  has  the  same  angular  momentum.  When  the 
gases  flow  radially  inward,  the  tangential  velocity  must  increase  to  conserve 
angular  momentum  This  leads  to  the  formation  of  a  void  core  of  zero  temper¬ 
ature  and  pressure  which  is  a  necessary  consequence  if  the  tangential  velocity 
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reaches  the  maximum  Isentroplc  speed. 

The  more  complex  grains  such  as  the  star,  the  wagon  wheel  and  the 
dendrite  add  considerable  difficulty.  The  mass  generated  at  various  radii 
from  the  axis  would  not  possess  a  unique  value  of  angular  momentum.  Also, 
there  could  exist  variations  In  the  burning  rate  due  to  stagnation  of  the 
faster  moving  gases  on  the  fingers  of  the  grain  and  pressure  variation  in 
the  radial  direction  balancing  the  acceleration  forces. 

When  viscous  forces  are  considered,  the  ideal  vortex  Is  upset.  This 
occurs  because  the  tangential  shear  stress  Increases  as  radius  decreases. 
At  the  centerline,  which  Is  a  free  boundary,  there  can  be  no  shear  moment. 
Therefore,  a  solid  body  vortex  is  established  close  to  the  axis  which  ex¬ 
tends  out  some  distance  into  the  chamber  and  merges  with  the  free  vortex. 


3,1,1  Theory  of  Burgers 

In  1940  and  1948,  Burgers  (11,  12)  discussed  an  exact  solution  to  the 
Navler-Stokes  equations  for  the  three-dimensional  vortex  flow  of  an  incom¬ 
pressible,  viscous  fluid. 

The  basis  for  this  solution  rests  upon  the  assumption  that  the  axial 

velocity  varies  with  axial  distance  only, 

w  s  a  »  z  (3,1) 


Then,  employing  the  continuity  equation, 


m  +  L  my, 

3Z  y 


(3,2) 


the  radial  velocity  becomes. 
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u  =  -  ar  +  K/r 


C3.3) 


where  K  is  a  constant  of  Integration  which  may  be  a  function  of  z  only. 

K  Is  then  eliminated  on  the  grounds  of  the  finiteness  of  the  radial 
velocity  at  the  centerline. 

u  =  -  ar  (3.4) 


The  radial  and  axial  momentun  equations  are  satisfied  by  the  above 
expressions.  The  tangential  momentun  equations  for  axlsymmetrlc  time 
independent  flow  is, 


2 

9 V  *«/  _  /  3 U  .  1  3 U 

-  Ar-  -  AV  -  + 

ar 

where  v  =  kinematic  viscosity. 


A  = 


vL 


Area 


=  const. 


(3.5) 


m/L  =  mass  per  unit  length 


The  following  solution  was  obtained  for  Eq.  (3.5)0 


V 


2rr  ^  ~  e 


■A  r  / 2r 


) 


r  =  2tt  V  •  r 


(3.6) 


where  r  is  the  circulation  about  the  axis. 
The  vortlcity  is  given  by 


(i) 


A  r  aAr2/2v 
2^e 


(3.7) 
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The  total  dissipation  per  unit  length  of  Z  is  given  by 


♦  = 


dr  2irry 


(3.8) 


It  can  be  seen  that  the  dissipation  is  independent  of  viscosity  and  pro¬ 
portional  to  the  third  power  of  the  tangential  velocities. 

Burgers  thus  obtained  an  exact  solution  for  a  three-dimensional 
vortex  flow.  The  results  are  completely  analytical  and  are  therefore 
quite  useful.  An  actual  flow  whose  axial  velocity  varies  linearly  with 
axial  position  while  the  tangential  and  radial  velocities  vary  only  with 
radial  position  may  not  be  physically  realizable.  However,  It  may 
approximate  a  psuedo  one-dimensional  flow  with  rotating  boundaries  through 
which  mass  is  Injected. 

Rott  (13,  14)  has  considered  Burger's  solution  in  more  detail  by 
allowing  the  velocity  profiles  to  vary  with  time.  He  also  has  calculated 
the  resultant  pressure  and  temperature  distributions. 


3.1.2  Theory  of  Einstein  and  LI 

In  1951,  Einstein  and  Li  (15)  studied  the  flow  in  an  emptying  con¬ 
tainer  filled  with  a  liquid.  The  model  for  the  analysis  was  essentially 
that  of  an  incompressible,  viscous,  steady,  axlsymmetric  flow.  It  was 
assumed  that  the  average  axial  velocity  u  was  small  while  the  variations 
in  the  axial  velocity  could  be  significant.  The  flow  field  was  divided 
into  two  portions.  A  core  having  the  same  radius  as  the  drain  opening 
was  assumed  in  which  the  vol unmetric  flow  rate  was  given  by 

Q  =  Q0  r2/ro2 


(3.9) 
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where, 

rQ  =  drain  radius 

Q0  =  The  total  volumetric  discharge  rate. 

In  the  remainder  of  the  flow  field,  from  the  drain  radius  to  the 
outer  walls, 


Q  =  Q0  =  const. 

R  =  The  cylinder  radius 

In  addition,  It  was  assumed  that, 


(3.10) 


2  2 

u  =  u(r)  u  <<  v 

2  2 

v  =  v(r)  w  <<  v 

w  =  0  ;  r  >  rQ 


Einstein  and  Li  found  that  In  the  core  region, 


=  (1  -  eA/2(r/ro)2) 


o  -  <rA/2) 


(3.11) 


where, 


A  =  (Q0/2*Yv) 

vQ  =  The  tangential  velocity  of  the  cylinder 
v  =  The  kinematic  viscosity 


In  the  outer  region  the  tangential  velocity  Is  given  by, 


V  - 


V  = 


v„  K 

(_0_0) 

\  r  ' 


(v_i_R) 

'  r  ' 


A(1 


(A  -  2)(1  -  eA/2(r/ro)  ) _ 

■  e‘A/2  •  (£)A_2)  -  2(1  -  e_A/2) 


-A/?  /’JL\  ~  ( A-2 ) 

A  -  2(1  -  e  A/2)(l  -  'R' 

A-2  5 ' 


A(1  -  e_A/2  £  (A-2))  -  2(1 


e-A/2) 


,  A  j*  2 


,  A  =  2 


(3.12) 


In  addition  to  these  laminar  flow  cases,  Einstein  and  Li  also  set  up 
the  equations  for  the  turbulent  motion  of  the  vortex.  They  concluded  that 
a  solution  equilivant  to  the  laminar  case  could  be  obtained  If  the  con¬ 
stant  A  Is  redefined  to  be, 


Ae  =  Q0/2tte(5  +  v) 

where, 

er|f(v/r)  =  -  (uV) 


(3.13) 


(3.14) 


3.1.3  Theory  of  Donaldson 

In  a  Ph.D.  thesis  (16)  and  other  reports  (17,  18),  Donaldson  presents 
a  rather  complete  review  of  the  solutions  for  the  Navier-Stokes  equations 
for  viscous,  incompressible,  driven  vorticies  which  have  been  accumulated 
over  the  years.  He  also  presents  a  great  deal  of  the  theory  and  equations 
for  such  motions  in  general  form. 

His  contribution  to  the  theory  of  viscous  vortex  motion  is  the  study 
of  the  Navier-Stokes  equations  for  a  three-dimensional  vortex  whose 
tangential  velocity  component  is  a  function  of  the  radial  coordinate  only. 

Under  the  above  catagory,  two  families  of  solutions  were  studied. 
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1 .  w  s  k  «  f(r) 

(3.15) 

2.  w  =  k  •  z  •  f(r) 

(3.16) 

The  following  physical  Interpretation  of  the  above  families  of  solutions 
were  given.  The  first  family  is  a  limited  one  since  any  mass  added  at 
the  boundaries  must  be  absorbed  by  a  sink  at  the  centerline.  However, 
the  second  family  where 

w  =  k  *  z  •  f(r)  (3.17) 


may  be  Interpreted  as  a  pseudo,  two-dimensional  problem  where  mass  Is 
added  through  a  rotating  container  and  Is  uniformly  distributed  In  the 
flow  field. 

When  the  axial  velocity  Is  of  the  form 

w  =  k-|  •  z  •  f.|(r)  +  k2  f2(r)  (3.18) 

Donaldson  shows  that  the  only  consistent  forms  of  the  tangential  and 
radial  velocities  are. 


u  =  u(r)  =  w  g  (x) 
v  =  v(r)  -  VQ  h  (x) 

where, 


(3.19) 


VQ  =  The  tangential  velocity  at  the  wall 


g(x)  =  An  arbitrary  function  of  radius  for  the  radial  velocity 
h(x)  =  An  arbitrary  function  of  radius  for  the  tangential  velocity 
R  =  The  radius  of  the  rotating  cylinder 
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Choosing  only  the  first  portion  of  the  general  expression  for  the 
axial  velocity, 

Xf(>0  (3.20) 

and  substituting  the  above  expression  Into  the  continuity  equation  gives, 

^  Hx  (xg)  =  -  f  •  (3.21 ) 

The  general  expression  for  the  velocities  may  then  be  substituted 
Into  the  axial  momentum  equation, 

fi--r(f)  t-9f'  -  +  +  ■  <3-22> 

Integrating, 

P  =  1/2  eJ  (f)2  [-  gf'  -  f2  +  ^  (f"  +  I')]  +  t(x)  (3.23) 

where, 

1 (x )  =  a  constant  of  partial  Integration 
Nw  =  WR/v  . 

Substituting  the  general  velocity  equations  into  the  radial  momentun 
equations  shows  that 

=  function  of  x  alone  .  (3.24) 

o  r 

Therefore, 

-  gf'  -  f2  +  ^  (f"  -  |-)  =  const.  =  c  .  (3.25) 

In  which  case. 
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|L  =  o  +  — =  function  of  x  alone  .  (3.26) 

dr  ar 

The  function  g(x)  may  be  eliminated  from  Eq.  (3.25)  by  substituting  the 
continuity  relationship,  Eq.  (3.21) 

f'f"'  _  (f"  _  f 7x)[f"  -  Nw(f2  -  c)]  -  Nw  f (f 1  )2  =  0  (3.27) 

Donaldson  sets  up  a  general  solution  for  the  differential  equations 
given  In  Eq.  (3.27).  For  the  case  of  an  Infinite,  rotating,  porous 
cylinder,  the  boundary  conditions  are, 

f(o)  =  1,  f 1 (o)  =  0,  f(l)  =  0  (3.28) 

In  which  case, 

C  =  -  (R/z)  •  (R/pu.2)  •  (f£)  . 

The  constant  C  may  be  determined  from  the  boundary  conditions  as  follows. 

C  =  f2(o)  -  RJ  +  IIS  •  (3.29) 

By  examining  the  limiting  process  In  Eq.  (3.29), 

(^^)  =  f"(o)  .  (3.30) 

the  following  result  was  obtained. 

C  =  1  -  2/Nw  f"(o)  .  (3.31) 

Now  the  Eq.  (3.27)  can  be  written, 

f'f"  -  (fn  -  f'/x)[f"  -  2f"(o)]  +  Nw(f2  -  l)(f"  -  f'/x) 

-  Nw  f(f')2  =  0 


(3.32) 
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Donaldson  Investigated  a  number  of  limiting  solutions. 


2.  N  =  +  *  2a,  Blowing  at  the  walls 

W 

3.  N  =  -  «  3a.  Suction  at  the  walls 

1 .  When  Nw  =  0,  v  =  « 

f'f"  -  (f"  -  f 7x)[f"  -  2f"(o)]  (3.33) 

It  was  found  that  there  are  an  infinite  number  of  solutions. 

f(x)  =  1  -  x2"1  m  =  1,  2,  3  ....  (3.34) 

However,  the  only  physically  possible  solutions  is, 

f(x)  =  1  -  x2  (3.35) 

2.  When  Nw  =  +  v  =  0 

(f2  -  l)(f"  -  f'/x)  -  f(f'2)  =  0  (3.36) 

In  this  case,  the  number  of  boundary  conditions  must  be  reduced, 
since  the  second  order  derivatives  were  eliminated, 

f (o )  =  1  ;  f'(o)  =  0  (3.37) 

The  general  form  of  the  solution  in  this  case  is 

f(x)  =  cos(ax2)  (3.38) 


Since  the  viscosity  has  been  eliminated,  the  axial  velocity  may  have  any 
desired  value  at  the  bounding  walls,  except  when  blowing  occurs,  in  which 
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case  the  axial  velocity  is  zero.  The  simplest  of  all  solutions  is  that 
when  a  =  1.  Then,  f  =  1  for  both  suction  and  blowing. 

2a.  Blowing  at  the  Walls  N  =  -  <=,  v  =  0 

_ _ W _ _ 

Since  In  blowing,  the  flow  Is  aware  of  the  presence  of  a  boundary, 

the  constant,  a,  must  be  chosen  so  that  the  axial  velocity  Is  zero  at  the 

wall  despite  the  vanishingly  small  viscosity.  Thus,  there  may  be  patterns 
of  alternating  forward  and  reverse  flow  for  the  axial  velocity. 

w  =  W  jj-  cos[(^l)  „x2]  (3.39) 

2b.  Suction  at  the  Walls  Nw  =  “  "» v  ~  0 

f  =  cos(ax2)  (3.40) 

In  this  case  the  constant,  a,  may  be  chosen  to  have  any  value  desired. 
In  general,  however,  the  form  of  the  solution  Is  the  same  as  for  the  blowing 
at  the  wall  case  except  that  the  requirement  for  the  axial  velocity  to 
vanish  at  the  wall  may  be  removed.  One  additional  point  arises,  however, 
because  when  the  magnitude  of  the  suction  is  large,  the  pressure  gradient 
may  be  adverse  to  the  net  mass  flow  so  that  any  number  of  Inward  and  out¬ 
ward  flow  cells  may  be  present.  Some  of  the  results  of  these  studies  are 
presented  In  Fig.  5. 


3.1.4  Theory  of  Lay 

A  paper  was  published  by  Lay  (19,  20)  In  1959  for  the  flow  In  a  tube 
with  superposition  of  the  free  vortex  onto  a  uniform  axial  flow.  The 
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analysts  Is  strictly  applicable  for  Inviscid,  irrotational ,  flows  since 
velocity  potentials  are  employed .  However,  compressibility  effects  are 
fncl uded* 

The  familiar  ax (symmetric  velocity  potential  and  stream  function 
equations  are  derived  for  the  case  of  tangential  and  radial  velocities 
only.  These  two  sets  of  equations  are  then  combined  under  the  hodograph 
transformation. 


V  =  q  COS0  s  II” 
u  s  q  sine  =  ||= 

The  equations  are  then.  In  terms  of  the  unknowns, 


(3.41) 


if/  -  Ip(q,e) 

where, 

q  -  the  velocity 

e  -  angle  of  the  velocity  vector 


(3.42) 


The  superposition  of  the  axial  flow  is  achieved  by  assuming  a  sink  flow 
in  two  dimensions  combined  with  a  uniform  axial  velocity  which  does  not 
alter  the  sink  flow  pattern.  It  was  possible  to  arrive  at  the  sink  flow 
vortex  flow  field  from  the  velocity  potential  and  the  stream  function 
equations  when  transformed  by  the  hodograph  transformation.  This  sink 
flow  vortex  was  then  uniformly  translated  by  superposition  of  a  constant 
axial  velocity. 

In  conclusion,  the  direction  of  the  streamlines  were  determined  for 
the  superposition  of  a  sink  flow  vortex  under  uniform  translation.  The 
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position  of  the  radius  of  maximum  isentropic  speed  was  also  calculated. 

3.1.5  Theory  of  Pengelly 

Pengelly  (21)  in  1957  presented  an  analysis  of  a  two-dimensional, 
compressible,  viscous  vortex.  The  radial  velocity  was  assuned  to  be  small 
compared  to  the  tangential  velocity  and  heat  transfer  was  neglected. 

The  model  for  the  flow  consists  of  a  rotating,  porous  container.  The 
gases  leave  the  cylinder  wall  with  a  combined  radial  and  tangential  velocity. 
A  schematic  of  the  physical  model  for  this  theory  is  presented  in  Fig.  6. 

The  gases  then  spiral  inward  to  a  radius  where  the  viscous  stresses  are 
zero,  at  which  point  the  radial  flow  is  converted  to  an  axial  flow  and 
exhausted  out  of  the  chamber.  Another  possibility  lies  in  the  removal 
of  the  Incoming  mass  by  slots  at  the  end  wall.  Thus,  the  flow  field  can 
be  treated  as  a  truly  two-dimensional  flow.  Therefore,  the  mass  flow  rate 
Is  given  by, 

m  =  puL2wr  (3.43) 

where 

L  =  the  length  of  the  chamber. 

All  gradients  with  respect  to  z  and  e  are  zero.  Therefore,  the 
following  form  of  the  Navier-Stokes  equations  can  be  written. 

a(du/dr  +  v/r)  +  rd2v/dr2  +  (dv/dr  -  v/r)  (3.44) 

dP/dr  =  (pv2/r)  {1  +  (u2/v2)(l  +  ^L)  (3.45) 


where 


FIGURE  6  ANG.  VELOCITY  WITH  RESPECT  TO  RADIAL  POSITION 
FOR  SEVERAL  VALUES  OF  FLOW  RATIO  ,o  . 

(Taken  from  reference  (21)  ,  Pengelly) 


In  addition,  the  energy  equation  Is  given  by 


dTs  .  1  rl  dP  r±.  dv  /,  ,  u2v  ,  v2  d  ,uv2i 

3 r“-  3^  V  3F-  ^  +  v  dF  (1  +  +  2  dF  (7J  ' 


(3.46) 


where, 


=  the  dlssapatlon  function 


*  ■  V#2  [f^-  n2  -  <^)2  {(^f)2  -  n  ]  (3.47) 


T$  =  the  stagnation  temperature. 


The  above  equations  are  exact  under  the  listed  assumptions  for  the 
case  of  a  two-dimensional  vortex.  Equation  (3.44)  can  be  Integrated 
without  approximation  for  the  case  of  o  equal  to  a  constant  to  give  the 
following  results. 

r  =  A  -  B/r°~2  a  t  2 
r  =  A2  -  B2  lnr  o  =  2 
v  =  r/2ur  . 


The  remaining  conservation  equations  can  be  integrated  provided  that 

2  2 
u  <<  v  , 

so  that, 

♦  -  »(*>*  {*$£  -  U*  (3.48) 

The  results  of  a  numerical  example  show  that,  under  certain  conditions, 
the  action  of  the  viscous  forces  can  decrease  the  stagnation  temperature 
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near  the  centerline  as  Is  found  in  the  Ranque-HIlsch  vortex  tube.  It  should 
be  emphasized  that  the  theory  is  applicable  only  to  a  two-dimensional  vortex 
with  laminar  flow.  The  mass  input  is  exhausted  out  of  the  cylinder  either 
through  slots  in  the  end  walls,  or  axially  through  some  small  hole  on  the 
centerl ine. 

Pengelly  also  discusses  the  possibility  of  applying  the  above  analysis 
to  the  turbulent  case.  He  concluded  that  it  may  be  possible  if  some  effective 
turbulent  viscosity  can  be  determined  which  can  be  employed  In  the  radial 
Reynolds  nunber. 


3,1.6  Theory  of  Deissler  and  Perlmuter 

In  1958,  Deissler  and  Perlmuter  (22)  presented  an  analysis  of  the  flow 
in  a  vortex  tube.  The  model  chosen  for  study  was  incompressible,  steady, 
viscous,  three-dimensional  flow.  It  was  assumed  that  the  tangential  and 
radial  velocities  were  functions  of  radius  only.  Then  it  was  shown  that 
only  a  linear  variation  of  axial  velocity  with  axial  position  is  consistent 
with  these  assumptions.  (See  Fig.  6). 

The  variation  in  axial  velocity  was  assumed  to  be  uniform  In  both  a 
core  region  and  an  outer  region.  However,  a  step  variation  in  the  velocities 
was  allowed  between  these  regions. 

Based  upon  this  model,  the  authors  were  able  to  determine  the  form  of 
the  tangential  velocity  as  a  function  of  radius  and  several  other  parameters. 
The  complete  tangential  momentum  equation  is  given  by. 


pu 


3_V 

ar 


3v  .  uv 

PUI  T —  +  p  — — 

M  3z  v  r 


3__  /  avt  a  r  ,av 
az  ^yazJ  ar  Ly  ^ar 


:)] 


+  i ]L  (11.1) 

r  var  r' 


(3.49) 


Under  the  assumption  that  V  =  v(r)  and  u  =  constant,  Eq.  (3.49)  can  be 
written  as. 


„  3v,  uv  _  ,d  v  ,  1  dv  V  . 

pU  W*  0  —  -  "  r  3F'  r*> 


(3.50) 


Therefore,  since  the  product  pu  Is  a  function  of  radius  only,  the  general 
solution  Is, 


V 


dr  +  C2/r 


(3.51) 


Where  C2  =  0  for  finiteness  of  velocity  at  r  =  0.  The  particular  solutions 
for  the  core  region  and  the  outer  region  are  then  determined  by  nunerlcal 
integration. 

In  addition  to  the  above  kinematic  study,  the  energy  equation  was 
also  solved  assuming  that  the  temperature  is  a  function  of  radius  only. 

The  solution  for  tangential  velocity  and  temperature  are  given  as  a  function 
of  radius  in  terms  of  the  following  parameters. 


r  /r,  =  the  ratio  of  the  assumed  core  radius  to  the  cylinder  radius 
wc/w  =  the  ratio  of  the  core  mass  flow  to  the  total  mass  flow 
Rg  =  the  radial  Reynolds  nunber 

The  results  indicate  that  with  high  Reynolds  numbers  the  flow  becomes 
increasingly  in  the  nature  of  a  free  vortex.  However,  there  always  remains 
a  small  region  near  the  centerline  of  solid-body  vortex  motion  which  grows 
as  the  Reynolds  number  decreases.  The  effect  of  wc/w  is  not  great  when 
the  Reynolds  number  is  large.  When  wc/w  is  small  and  the  Reynolds  nunber 
is  small,  the  flow  tends  to  be  more  free  vortex  in  nature. 
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The  temperature  distributions  show  that  there  is  very  little  energy 
separation  except  at  the  lowest  Reynolds  numbers.  This  suggests  that  the 
observed  energy  separation  in  the  Ranque-Hilsch  tube  experiments  must  be 
attributed  to  turbulence  effects  which  would  essentially  increase  the 
apparent  viscosity  of  the  flow.  Therefore,  a  turbulent  model  was  proposed 
based  upon  the  additional  assumption  of  a  constant  eddy  viscosity.  This 
made  no  change  in  the  velocity  distribution  for  a  given  Reynolds  number. 
For  the  calculation  of  the  turbulent  temperature  distribution,  it  was 
necessary  to  make  the  following  assumptions. 

p  -  pe 

k  1  p  cp 

e  ”  eH 

An  additional  term  was  then  added  to  account  for  the  effect  of  the  radial 
pressure  gradient  upon  the  turbulent  heat  transfer. 

In  conclusion,  the  turbulent  temperature  analysis  was  shown  to  agree 
with  experimental  vortex  tube  data  when  the  proper  choice  of  the  eddy  vis¬ 
cosity  was  made.  Essentially,  the  effects  of  the  Increased  viscosity  and 
the  radial  pressure  gradient  upon  the  turbulent  heat  transfer  allowed  a 
greater  energy  separation,  thus,  a  better  correlation  with  experimental 
measurements. 


3.1,7  Theory  of  Sullivan 

In  1959,  Sullivan  (23)  reported  upon  an  analog  to  Burger's  solution 
(ref.  11,  12)  which  is  characterized  by  the  exl stance  of  two  cells.  In 
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this  model,  the  flow  does  not  spiral  toward  the  axis  then  out  along  It, 
but  has  a  central  region  of  reverse  flow.  This  solution  came  about  during 
the  course  of  an  Investigation  of  the  complete  family  of  solutions  based 
upon  velocity  profiles  of  the  following  form  in  an  Incompressible,  viscous, 
steady,  axlsymmetrlc  flow  field. 

u  =  u(r) 

v  =  v(r)  (3.52) 

w  =  k  •  z  •  w(r)  , 


The  particular  solution  In  this  report  Is  for  the  case  of  an  Infinite 
rotating  porous  container  through  which  liquid  Is  flowing.  The  radial 
and  tangential  velocities  and  their  derivatives  are  required  to  vanish 
at  the  centerline.  At  the  container  walls  the  velocity  of  the  boundary 
and  the  tangential  velocity  are  the  same. 

The  two  cell  solution  solution  arises  from  the  case  where  the  radial 
Reynolds  number  approaches  minus  Infinity.  Physically,  this  means  that 
large  Inward  radial  velocities  exist  at  the  container  walls  so  that  the 
inertial  forces  are  much  larger  than  the  viscous  forces.  The  following 
solution  was  obtained. 


2 

u  =  -  ar  +  6  (v/r)[l  -  e"(ar  /2v)] 
w  =  2az  [1  -  3  e"(ar2/2v)] 


v  =  r/2irr  [H(arZ/2\j)/H(«)] 

H(x)  =  X  exp  {-  5  +  3  I*  dT}dt 


(3.53) 

(3.54) 

(3.55) 

(3.56) 
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THE  TANGENTIAL  VELOCITY  PROFILES 
BURGERS’ SOLUTION  SULLIVAN^  SOLUTION 


FIGURE  7  A  COMPARISON  OF  THE  THEORY  OF  BURGERS  AND  SULLIVAN 
(Taken  from  reference  (27),  Lewellen) 
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Sullivan  compares  this  solution  to  that  of  Burger  and  finds  that, 
although  the  velocity  expressions  are  different,  the  axial  pressure 
gradients  and  dissipation  are  the  same.  Thus, 


3P/3t  =  -  4  a  pz 

,2 


(3.57) 


$  ■ 


_  p  A  r‘ 


4  TT 


The  radial  pressure  function  is,  In  general,  quite  different, 

P  =  PQ  -  p/2  {4  a2z2  +  a2r2  +  36(\>2/r2)[l  -  ear2/2v]2}  + 

+  p  f  (v2/r2)dr 
■^0 

whereas  Burgers  and  Rott  (11,  12,  13,  14)  found, 

rr 


(3.58) 


P  -  P0  -  (p/2)  [4  a2z2  +  a V]  + 


(u/r)  dr  . 


(3.59) 


3.1.8  Theory  of  Slbulkln 

In  a  series  of  three  articles  (24,  25,  26),  Slbulkln  presents  a  method 
of  analyzing  certain  vortex  flows.  The  basis  of  the  method  Is  to  establish 
the  unsteady  equations  of  tangential  and  radial  momentun  for  viscous  flow 
and  the  energy  equation  allowing  radial  transfer  of  energy  only. 

Radial  Momentun 

pro)2  =  3P/3r  (3.60) 

Tangential  Momentum 

3 o>  _  \>  3 _  /  3  3u)\ 

3 1  "  3  3r  3r; 
r 


(3.61) 
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Energy 

dT 

T 


The  restrictions  on  these  equations  are: 

1.  M  «  1 

2.  u  -  0 

3.  Only  radial  energy  transfer  Is  allowed 

4.  Laminar  flow 

5.  Axial  and  radial  shear  forces  neglected. 


=  c. 


3t  at 


_  1  3P 


±iUL+  v  a_  rr 
p  3z  r  3r  L 


3  i ill: 

ar  - 


pr  3r  v  dr' 


(3.62) 


The  Initial  conditions  for  the  problem  are  that  of  given  Initial  property 
profiles  determined  from  experiment. 

The  analytical  approximation  to  these  Initial  conditions  are  a 
quiescent  core  of  zero  temperature,  pressure,  and  velocity,  with  a  free 
vortex  distribution  around  It.  These  Initial  conditions  are  gradually 
smoothed  out  as  the  Integration  In  time  proceeds.  The  boundary  conditions 
In  time  and  radius  are: 


f-  (0,t)  =  0 

u  (R,t)  =  0 


3 1 

(0,t)  =  0 
3i. 

3^  (R.t)  =  0 


(3.63) 


t ^  *  the  total  temperature. 


These  boundary  conditions  are  for  a  stationary  tube  which  extends 
to  Infinity.  The  axial  position  then  is  some  function  of  time.  These 
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equations  were  solved  simultaneously  by  finite  difference  methods.  The 
results  indicate  that  the  tangential  velocity  decays  in  time  while  the 
pressure  gradient  of  the  imposed  inlet  profile  Is  gradually  smoothed. 

The  total  temperature  peak,  however.  Indicated  a  gradual  tendency  to 
move  inward.  The  initial  conditions  and  some  results  are  presented  in 
Fig.  7. 

3.1.9  Theory  of  Lewellen 

In  1962  and  1964,  Lewellen  (27,  28)  investigated  the  properties  of 
real  vortices  contained  In  cylinders  with  small  axial  flow.  The  type 
of  flow  considered  is  one  with  strong  circulation  wherein  the  radial 
velocities  are  small  compared  to  the  tangential  velocity  as  Illustrated 
in  Fig.  8.  These  solutions  are  compared  to  those  of  Donaldson  (16), 
Sullivan  (23),  Burgers  (12),  and  Rott  (13).  The  above  authors  Investigated 
exact  solution  to  the  Navler-Stokes  equations  when  the  following  conditions 
hold. 


W(r,z)  =  K 1  *  z  =  f-j(r)  +  K2  f2(r)  (3.64) 

where, 

u  =  u(r) 
v  =  v(r)  . 

These  analyses  may  be  contrasted  to  those  of  Einstein  and  LI  (15),  and 
Delssler  and  Perlmutter  (22),  where  approximate  techniques  were  employed. 

The  assumptions  leading  to  solutions  of  the  three-dimensional  flow 
in  the  analyses  listed  above  are  that, 
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2  2 

w  <<  v 

2  2 

u  <<  v 

(3.65) 

u  -  u( r) 
v  =  v(r)  o 


These  approximations  are  justified  when  the  axial  velocity  of  the  flow 
is  small  and  the  circulation  is  large.  In  the  present  paper,  Lewellen 
attempted  to  bridge  the  gap  of  these  two  methods  of  attack  and  to  point 
out  the  consequences  of  each  approach.  He  took  the  complete  Navier- 
Stokes  equations  for  axi symmetric,  viscous,  laminar,  incompressible, 
steady  flow  in  cylindrical  coordinates  and  employs  the  concept  of  the 
stream  function  to  eliminate  the  continuity  equation.  The  remaining 
equations  were  then  written  in  terms  of, 


r  =  vr  =  r(r,z) 

ip  -  ip(v,z)  =  The  stream  function 


(3.66) 


where  pressure  is  eliminated  by  cross  differentiation  of  the  axial  and 
radial  momentun  equations,  Lewellen  then  derived  the  equations  pertinent 
to  the  analyses  of  Donaldson  and  Sullivan  (16,  23).  He  showed  that  the 
solutions  given  by  Donaldson  and  Sullivan  are  not  applicable  to  the 
analysis  of  the  radial  sink  flow  to  a  small  opening  in  a  cylindrical 
tank.  However,  by  assuming  pertibations  in  the  Donaldson  and  Sullivan 
analysis  to  include  small  changes  of  ^  with  axial  position,  Lewellen 
showed  that  when  the  circulation  is  large,  there  can  be  large  changes  in 
the  stream  function  ^ . 
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To  Investigate  this  case,  Lewellen  expanded  the  expressions  for 
the  circulation  and  the  stream  function  In  the  following  manner. 

r  s  l  r  (n,c)en  where  n  =  (r/r  )2 
nso 

(3.67 

*  =  l  (n.c)en  z  =  z/l 

n=o 

The  solutions  thus  obtained  are  valid  for  variations  of  the  radial 
pressure  gradient  and  axial  velocities  and  therefore  are  not  necessarily 
specified  a  priori.  However,  the  solutions  are  quite  dependent  upon  the 
boundary  conditions  of  the  stream  function  which  In  physical  cases,  must 
be  related  to  the  boundary  layer  phenomenon  associated  with  the  end  section 
of  the  container.  It  is  noted  that  in  most  real  cases,  the  boundary  layer 
must  be  solved  to  give  proper  initial  conditions  for  the  calculations,  and 
some  consideration  should  be  given  to  the  recirculation  induced  by  the 
exhaust. 


3.1.10  Lewellen  1965 

In  this  paper,  as  In  earlier  papers,  Lewellen  (29)  investigated  the 
flow  associated  with  uniform  tangential  Injection  with  strong  circulation 
and  a  sink  on  the  axis  with  a  radius  smaller  than  that  of  the  container. 

The  theory  Is  developed  by  considering  a  linearized  solution  about  weak 
circulation  and  another  about  strong  circulation. 

A  particular  example  Is  discussed  in  which  the  swirling  flow  is 
generated  by  a  rotating  porous  container.  The  linearized  theory  for 
strong  circulation  with  weak  axial  flow  is  used.  It  Is  also  assumed  that 
the  Reynolds  number  associated  with  rotation  is  so  large  that  the  principal 
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variation  of  the  axial  velocity  occurs  In  the  end  wall  of  the  boundary 
layer,  while  the  radial  variation  occurs  In  the  cylinder  wall  boundary 
layer.  Lewellen  showed  the  large  swirl  forces  of  the  flow  to  be  essentially 
two-dimensional  except  In  a  region  near  the  boundary  layers  where  the 
principle  variation  In  the  velocities  occur. 


50 


4a  VISCOUS  BOUNDARY  LAYERS  IN  VORTEX  FLOWS 
4*1  Introduction 

In  the  previous  section,  theories  were  presented  which  may  be 
applicable  to  the  flow  inside  the  chamber  of  an  Internal -burning  rocket 
motor,.  Even  though  viscous  effects  were  in  general  included,  the 
Interaction  of  the  primary  flow  with  the  end  wall  of  the  container  were 
not  considered*  The  purpose  of  this  section  is  to  present  some  of  the 
considerations  pertinent  to  viscous  boundary  layers  around  rotating 
disks  and  In  nozzles*  The  bulk  of  this  theory  is  applicable  for  Incom¬ 
pressible,  laminar  flow  only*  However,  some  considerations  are  given 
to  turbulent  effects* 

4*2  Boundary  Layers  in  Nozzles 

The  analyses  presented  in  this  section  are  characterized  by  an 
Irrotational  core  of  an  Incompressible  fluid  which  Is  surrounded  by  a 
boundary  layer  formed  at  the  stationary  walls  of  a  conical  nozzle*  In 
1940,  Taylor  (30)  presented  an  analysis  of  the  swirl  atomizer*  This 
device  had  been  studied  previously  by  many  authors*  However,  Taylor 
Introduced  the  concept  of  an  irrotational  core  surrounded  by  a  viscous 
layer  attached  to  the  stationary  nozzle  walls*  Based  upon  measurements 
of  the  discharge  coefficients  of  swirling  nozzle  flows,  he  concluded 
that  the  boundary  layer  could  occupy  a  significant  portion  of  the  nozzle 
flow  area*  Thus  he  concluded  that  the  treatment  of  the  swirl  atomlxer 
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by  perfect  fluid  theory  was  not  accurate.  A  schematic  of  the  swirl 
atomizer  and  the  model  for  the  analysis  of  the  boundary  layer  In  nozzles 
with  swirling  flow  is  presented  In  Fig.  9. 

4.2.1  Theory  of  Taylor 

Later  In  1945,  Taylor  (31)  presented  a  more  complete  analysis  of 
the  boundary  layer  phenonema  In  the  nozzle  of  a  swirl  atomizer.  The 
following  list  of  assumptions  are  common  to  this  analysis  and  succeeding 
analyses. 

1.  The  boundary  layer  Is  thin. 

2.  The  pressure  variation  normal  to  the  boundary  layer  Is  small, 

3.  There  Is  an  Irrotational  core  which  imposes  Its  pressure  on 
the  boundary  layer. 

4.  Terms  Involving  radial  velocity  squared  are  neglected  with  respect 
to  tangential  velocity  squared. 

The  problem  Taylor  chose  to  solve  was  that  where  the  axial  velocity 
of  the  core  could  be  considered  negligible.  Taylor  used  the  following 
velocity  profiles. 


Rsma  5 

f(n)  -  (n  -  2n^  +  n^) 

(4.1) 

Rsma  * 

<f>(n)  =  (2n  -  n^) 

(4.2) 

.  R(a-8) 

6 

(4.3) 

6  =  The  boundary  layer  thickness 
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FIGURE  9  THE  BOUNDARY  LAYER  IN  A  CONVERGING  NOZZLE  IN 
SPHERICAL  POLAR  COORDINATES 
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a  =  The  cone  half  angle 
0  =  The  flow  angle. 

Taylor  derived  the  equations  for  Incompressible,  viscous  flow  In  the 
spherical  polar  coordinate  system  subject  to  the  general  assumptions  listed 
above.  The  equations  were  placed  In  Integral  form.  The  following  boundary 
conditions  were  satisfied  In  the  intergal  formulations. 

f-|4>=0;  n  =  0:  f  =  0,  4>  =  1 .  n  =  1  (4.4) 

Taylor  was  then  able  to  determine  the  growth  of  the  boundary  layer 
and  the  direction  of  the  flow  at  the  nozzle  surface.  In  conclusion,  two 
examples  were  presented  for  a  cone  half-angle  of  45°.  When  the  Injection 
velocity  was  45  meters  per  second  and  the  chamber  pressure  was  ten  atmospheres, 
the  boundary  layer  thickness  was  11  per  cent  of  the  total  radius  of  the  throat. 
When  the  chamber  pressure  was  one  atmosphere,  the  boundary  layer  was  fully 
developed  and  occupies  the  entire  flow  area.  A  plot  of  the  variation  of  the 
boundary  layer  thickness  with  respect  to  the  reduced  nozzle  radius,  R/RQ, 
is  presented  in  Fig,  10, 

4.2.2  Theory  of  Binnle  and  Harris 

In  1950,  Binnle  and  Harris  (32)  published  an  analysis  of  the  swirl 
atomizer.  In  this  paper,  some  of  the  general  theory  of  Irrotational ,  In¬ 
compressible  vortices  was  discussed.  They  also  presented  an  analysis  of 
the  boundary  layer  growth  in  the  nozzle.  The  assumptions  employed  were 
essentially  those  used  by  Taylor  except  that  the  axial  velocity  In  the 
core  was  not  necessarily  considered  to  be  neglegibly  small. 


(a)  SWIRL  ONLY;  (b)  STREAM  IN  (i  ONL?;  (c)  SWIRL  AND  STREAMING 


RESULTS  OF  A  SAMPLE  CALCULATION  BY  BlNNIE  AND  HARRIS  (32) 


FIGURE  10  SOME  RESULTS  FOR  THE  BOUNDARY  THEORY  IN  CONICAL  NOZZLES 
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By  assuming  the  following  velocity  profiles,  an  estimate  was  obtained 
for  the  boundary  layer  thickness,  6,  as  a  function  of  the  distance  In  the 
nozzle. 


w  = 

G(R)  4> ( n )  $  = 

2n  -  n2 

(4.5) 

V  = 

R  sin.  +(n)  "  3 

R(a-e) 

6 

(4.6) 

6  = 

The  boundary  layer  thickness 

(4.7) 

a  - 

The  cone  half  angle 

e  = 

The  flow  angle 

In  conclusion,  Blnnle  and  Harris  found  two  simultaneous  differential 
equations  Involving  the  pressure  and  velocity  profiles.  These  equations 
were  Integrated  using  a  forward  marching  numerical  method  and  assuming  con- 
tant  coefficients  over  small  steps.  The  results  of  their  calculations  were 
compared  to  Taylor's  analysis  (31)  which  was  based  upon  zero  actual  velocity 
In  the  core.  Binnie  and  Harris  found  that  the  result  for  the  boundary  layer 
thickness  Is  of  the  order  of  one-third  to  one-fifth  of  Taylors'  boundary 
layer  thickness  as  shown  in  Fig.  10. 

4,2=3  Theory  of  Cooke 

In  1951,  Cooke  (33)  published  an  article  concerning  the  boundary  layer 
buildup  of  Taylor's  problem  (31)  wherein  the  axial  velocity  of  the  core  was 
assigned  to  be  negligible.  The  main  purpose  of  this  article  was  to  show 
that  the  boundary  layer  thickness  associated  with  the  axial  velocity  was 
not  necessarily  the  same  as  the  boundary  layer  thinkness  associated  with 
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the  swirling  component.  In  other  words,  he  postulated  that  two  boundary 
layer  thicknesses  are  required. 


W/WQ  =  .995  -  defines  the  location  of  5 


v/v0  *  -995  -  the  °f  4 


(4.8) 


Cooke  incorporated  the  boundary  layer  conditions  used  by  Taylor, 
allowing  an  additional  set  of  boundary  conditions  for  a. 


W  =  0  V  =  0 
2 

^Mr  =  r/v  R  sina  at  e  *  a  (4.9) 

30* 

3*V  _  0 

— 5-  -  U 

30 

W  =  0,  3aj/30  *  0  at  0  =  a  -  6/R 

V  =  TKT  *  If  =  0  at  0  =  a  -  A/R 
sina  30 

The  solution  for  the  boundary  layer  thicknesses  then  preceeded  in  the 
manner  of  Pohlhausen  by  assuming  polynomials  for  the  velocity  profiles 
through  the  boundary  layer.  Cooke  found  that  the  boundary  layer  thickness 
determined  In  this  manner  was  approximately  one-half  that  of  Taylor's 
problem,  and  therefore,  between  the  analyses  of  Taylor  (31)  and  Blnnle  and 
Harris  (32),  as  shown  In  Fig.  10. 


4.2.4  Theory  of  Weber 

In  1956,  Weber  (34)  presented  an  analysis  of  the  boundary  layer  In  a 
conical  nozzle.  Once  again,  the  boundary  conditions  of  Taylor's  problem 
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were  assumed.  However,  in  this  paper,  the  Pohlhausen  intergal  technique 
was  used  for  both  laminar  and  turbulent  flows.  For  the  laminar  case,  the 
results  are  identical  to  Taylor's  results  because  the  same  velocity  pro¬ 
files  and  boundary  conditions  were  assumed. 

In  the  turbulent  case,  the  one-seventh  power  law  for  the  turbulent 
boundary  layer  was  used  for  the  tangential  velocity  distribution. 

V  =  V6  n1/7  (4.10) 

Then,  due  to  a  lack  of  data,  Weber  assuned  the  axial  velocity  to  be 

W  =  Wfi  E  (n1/7  -  n)  (4.11) 

so  that  the  boundary  conditions  of  zero  velocities  could  be  satisfied  at 
the  wall.  The  assumption  of  this  type  of  velocity  profile  makes  It  Impossible 
to  determine  the  wall  stresses  since  the  gradients  become  Infinite  at  the 
wall.  Therefore  an  experimentally  obtained  correlation  was  used  for  the 
wall  shear  stress. 

to)  =  0.0225  pV7  (v/y)1/4  (4.12) 

The  stress  components  at  the  wall  were  determined  next  and  substituted 
together  with  the  assumed  velocity  profiles  into  the  momentum  Integral 
equations.  The  resulting  boundary  layer  thickness  was  found  to  be  approxi¬ 
mately  one-fourth  that  of  Taylor's  laminar  case,  as  shown  in  Fig.  10. 

4.2.5  Conclusions 

In  conclusion,  it  appears  that  the  theory  for  the  development  of  the 
boundary  layer  in  a  simple  conical  nozzle  for  incompressible  flow  is  not 
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well  understood ,  However,  the  boundary  layer  does  seem  to  occupy  a  portion 
of  the  nozzle  which  can  be  large  when  the  chamber  pressures  are  small  and 
the  axial  velocities  are  small ,  In  a  rotating  rocket  motor,  the  walls  of 
the  motor  would  be  rotating  in  direct  proportion  to  their  radial  position,. 

The  effect  of  moving  walls  would  be  to  reduce  the  velocity  difference  be¬ 
tween  the  core  and  the  walls,  thereby  reducing  the  boundary  layer  buildup, 

4,3  Boundary  Layers  Near  Rotating  Disks 

This  section  deals  with  a  class  of  problems  associated  with  the  flow 
of  a  viscous  fluid  bounded  by  rotating  disks,  T,  Von  Kamian  (35)  first 
attacked  this  problem  In  1921,  He  obtained  a  solution  for  a  steady,  viscous 
flow  In  a  semi-infinite  region  bounded  by  a  rotating  place.  The  solution 
was  exact  In  that  no  assumed  velocity  profiles  were  required.  However, 
the  resulting  differential  equations  required  numerical  solution, 

4,3,1  Theory  of  Batchelor 

In  1950,  Batchelor  (36)  extended  Von  Karman's  analysis  to  include  a 
flow  bounded  by  two  Infinite  rotating  disks.  In  general,  two  classes  of 
problems  were  studied?  one,  in  which  the  plates  rotated  In  the  same  direction, 
and  the  second,  in  which  the  plates  rotated  in  the  opposite  sense,  Batchelor 
worked  out  these  solutions  qualitatively  to  illustrate  the  pertinent  features 
of  the  flow, 

Batchelor  began  with  the  equations  of  steady,  viscous,  incompressible, 
axlsymmetric  flow.  The  only  restriction  placed  upon  the  equations  was  the 
assumption  that  the  velocity  of  the  fluid  flowing  normal  to  the  disk  was  a 
function  of  axial  position  only. 
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W  =  W(z)  (4.13) 

The  boundary  conditions  are  those  of  zero  velocities  relative  to  the 
rotating  plates 


w  = 

u  =  0  ; 

v  =  wr  at  z  =  0 

w  = 

u  =  0  ; 

v  =  r  at  z  =  d 

where 

(o  =  angular  velocity  of  the  first  disk 
y  =  angular  velocity  of  the  second  disk 

or,  y.|=  angular  velocity  of  the  fluid  far  from  the  first  disk 

Solutions  were  obtained  for  one  and  two  disk  systems  as  a  function 
of  y-i/u.  The  results  of  these  studies  are  presented  In  Fig.  11. 

4.3.2  Theory  of  Rodgers  and  Lance 

In  1949,  Rodgers  and  Lance  (37)  presented  additional  studies  of  a 
rotatlonally  symmetric  flow  of  a  viscous  fluid  In  the  presence  of  an  In¬ 
finite  rotating  disk.  In  their  analysis,  Rodgers  and  Lance  treated  the 
following  cases  employing  similarity  profiles. 

1.  The  fluid  and  the  disk  rotating  in  the  same  sense 

2.  The  fluid  and  the  disk  rotating  In  the  opposite  sense 

3.  The  fluid  rotating  as  a  solid  body  at  infinity. 

The  particular  property  profile  functions  employed  are  given  below. 
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u  =  r  n  F(c) 
v  =  r  n  G(jj) 

w  =  (vn)1/2  H(c)  (4.15) 

P/c  =  vfi  P0(C)  +  1/2  k  n2  r2 

fi  s  The  angular  velocity  of  the  disk 

These  equations  were  then  introduced  into  the  Navler-Stokes  equations 
yielding  a  set  of  four  ordinary  differential  equations. 

Physically  acceptable  solutions  to  these  equations  were  found  for 
all  values  of  the  ratio  of  fluid  angular  velocity  at  Infinity  to  the 
disk  angular  velocity  for  case  1  described  above.  In  each  Instance, 
the  boundary  layer  was  attached  to  the  disk  and  the  motion  approached 
solid  body  rotation  In  the  fluid  far  from  the  disk. 

In  case  2,  where  the  fluid  Is  rotating  In  the  opposite  sense  as  the 
disk,  there  existed  regions  of  unsteady  solutions  only.  The  only  correction 
for  such  cases  was  the  application  of  suction  to  prevent  the  boundary  layer 
leaving  the  disk.  It  was  concluded  that  when  the  flow  along  the  disk  was 
radially  inward,  the  slmillarlty  solution  employed  may  not  be  valid,  espe¬ 
cially  if  separation  is  likely  to  occur. 

4.3.3  Theory  of  Mack 

In  a  series  of  reports  by  Mach  (38,  39),  the  problem  of  the  Interaction 
of  a  primary  vortex  field  with  a  stationary  disk  was  studied.  It  had  been 
shown  previously  by  experiment  that  the  presence  of  the  end  walls  cause  a 
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radial  In-flow  due  to  the  Imbalanced  radial  pressure  gradient  in  the  boundary 
layer.  It  was  the  purpose  of  Mack's  reports  to  Investigate  methods  of  an¬ 
alytically  predicting  the  amount  of  in-flow  and  its  effect  upon  the  primary 
stream  in  the  chamber. 

Various  techniques  were  employed  to  obtain  solutions  for  the  incompress¬ 
ible,  boundary  layer  equations  on  a  disk,  The  method  of  Von  KarmSn  was  em¬ 
ployed  initially,  however,  Taylors'  (31)  and  Cooke's  (33)  methods  were  also 
applied.  Stewartson  velocity  profiles  were  used  In  these  integral  techniques 
and  these  results  were  compared  to  a  set  of  eleven  pairs  of  velocity  profiles. 
The  Stewartson  profiles  are  very  useful  In  this  respect  since  they  allow 
exact  solutions  of  the  boundary  layer  equations  at  the  edge  of  the  disk. 

To  complete  the  study  primary  flow  field  tangential  velocities  of  the 
following  form  were  assuned. 

V  a  l/rn 
-  1  <  n  <  1 

These  were  employed  to  obtain  certain  boundary  conditions  for  the  boundary 

layer.  In  the  case  of  n  =  +  1 ,  or  a  free  vortex,  the  radial  mass  flow 

Increase  monotomlcally  with  decreasing  radius.  For  n  t  1,  the  radial  mass 

flow  reaches  a  maximum  at  some  radius  less  than  the  centerline.  When  n  =  -1 , 

a  solid  body  vortex,  the  radial  mass  flow  is  only  30%  of  the  free  vortex 

case  and  Its  maximun  value  occurs  at  r/r  =  .74. 

o 

4.3.4  Theory  of  King  and  Lewellen 

King  and  Lewellen  (40,  41,  42)  have  prepared  a  number  of  reports  on 
the  action  of  the  viscous  boundary  layer  associated  with  the  end  walls  of 
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a  vessel  containing  a  vortex  flow.  These  studies  were  made  because  It  had 
been  found  experimentally  that  these  end-wall  boundary  layers  significantly 
affect  the  flow  In  the  chamber.  Previous  studies  of  vortices  In  chambers 
are  strlckly  applicable  for  Infinite  cylinders  only  since  boundary  layers 
on  the  end  walls  were  neglected.  It  should  be  noted  that  these  solutions 
are  limited  to  stationary  end  walls. 

The  calculations  were  performed  nunerlcally  employing  several  variations 
of  the  classical  Kanman-  Pohl hausen  Integral  technique  for  Incompressible 
flows.  The  tangential  velocity  of  the  main  stream  was  varied  by  assunlng, 

V  o  r^ 

where, 

-  1  <  rt  <  +  1 

The  value  of  the  exponent  n  was  varied  parametrically.  It  was  found 
that  the  most  drastic  effects  upon  the  main  stream  occurred  when  n  =  -1. 

Since  the  end  wall  was  a  stationary  boundary,  the  tangential  velocity  of 
the  free  stream  Is  retarded  as  It  approaches  the  wall,  thus  losing  cen¬ 
trifugal  force.  When  n  =  -  1,  the  primary  stream  was  a  free  vortex,  and 
the  pressure  gradient  associated  with  a  free  vortex  Is  Imposed  upon  the 
wall.  The  result  was  an  unbalanced  pressure  gradient  at  the  wall.  Thus, 
an  Inward  flowing  radial  flow  was  generated,  associated  with  the  unbalanced 
radial  pressure  gradient.  Continuity  required  the  presence  of  an  axial 
velocity  In  the  presence  of  the  Induced  radial  velocity.  Thus,  secondary 
flows  are  set  up  In  the  main  chamber.  In  addition,  the  boundary  layer  was 
able  to  control  a  great  deal  of  the  total  mass  flow  rate. 
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The  results  of  the  studies  performed  are  compared  to  the  theories  of 
Cooke  (33),  Taylor  (31),  Weber  (34),  and  Mack  (38),  It  is  found  that  when 
n  f  -  1 ,  the  results  of  all  the  analyses  were  approximately  the  same*  How¬ 
ever,  when  n  -  -  1,  considerable  discrepancies  in  the  values  of  the  radial 
inflow  were  found, 

4.3,5  Theory  of  Rosenweig,  Lewellen,  and  Ross 

In  a  paper  published  in  late  1964,  Rosenweig,  Lewellen,  and  Ross  (43) 
presented  an  analysis  of  a  contained  vortex  in  which  fluid  is  Injected 
tangentially  through  porous,  rotating  walls.  This  paper  differs  from  those 
published  previously  by  Lewellen  (27,  28)  in  that  the  flow  field  Is  divided 
into  three  fundamental  regions,  A  schematic  of  the  physical  model  is  pre¬ 
sented  In  Fig,  12. 

Region  I  consists  of  the  primary  flow  region  which  is  free  from  edge 
effects  and  boundary  layer  interaction.  It  occupies  a  volume  defined  by 
the  tank  walls,  the  end-wall  boundary  layer,  and  the  drain  opening.  Region 

II  consists  of  the  boundary  layer  associated  with  the  end  walls.  Region 

III  is  bounded  by  an  Imaginary  cylinder  extending  up  from  the  drain  opening. 
Each  of  these  regions  is  considered  to  be  viscous  in  nature.  However,  in 
Region  II,  axial  gradients  predominate,  whereas  in  Region  II  and  III,  the 
radial  gradients  predominate. 

Approximate  solutions  are  employed  In  each  of  these  regions.  A 
complete  flow  field  is  built  up  by  prescribing  suitable  matching  conditions 
between  these  regions.  The  solution  is  based  upon  the  assumption  of  large 
circulation  and  small  axial  velocities. 
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The  approximate  boundary  layer  solution  of  Rott  (44)  Is  used  along 
the  end-walls  and  Is  matched  to  the  solution  of  the  primary  flow  field. 

The  authors  specify  the  stream  function  distribution  In  the  central  core, 
region  III.  A  provision  is  made,  however,  to  allow  for  the  Injection  of 
fluid  from  the  end-wall  boundary  layer  Into  the  drain  opening.  The  equations 
apply  to  an  incompressible,  laminar  flow  field  where  appropriate  boundary 
layer  assumptions  are  made  in  Regions  II  and  III.  The  solutions  are 
parameterized  In  terms  of 

Rgr  =  the  radial  Reynolds  number 
re/rQ  =  the  ratio  of  drain  opening  to  the  cylinder  radius 

"t 9  the  dimensionless  axial  coordinate 
*  2 6  ro  ro 

A  =  — p  =  the  end-wall  boundary  layer  interaction  term. 

(Ret)  0 

where, 

rQ  =  the  cylinder  radius 


Tl  =  the  length  of  the  cylinder 
Q0  =  the  volume  flow  rate 
Ret  =  tangential  Reynolds  nunber 

Solutions  are  presented  in  terms  of  the  ratio,  r/rQ,  the  ratio  of  the 
point  circulation  to  the  wall  circulation,  as  a  function  of  the  reduced 
radial  coordinate,  r/rQ.  The  solutions  are  computed  for  constant  values 
of  the  above  parameters.  However,  only  rg/r0  =  1/6  Is  used.  It  Is  found 
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that  when  A  <  1»  the  boundary  layer  Interaction  Is  not  large.  The  above  Is 
also  true  when  Rgr,  Is  small,  regardless  of  the  magnitude  of  A.  The  axial 
variation  of  T/r  with  5  Is  small  In  these  cases.  However,  when  A  and,  or, 

R  are  large,  the  boundary  layer  Interaction  is  felt  throughout  the  entire 
flow  field.  Also,  the  effects  of  mass  ejection  from  the  boundary  layer  tends 
to  Increase  until  reverse  flow  and  recirculation  are  necessary  to  keep  the 
boundary  layer  supplied  with  fluid. 

4.3.6  Theory  of  Rott  and  Lewellen 

A  comprehensive  review  of  the  phenomena  of  boundary  layer  interaction 
with  a  vortex  flow  field  Is  presented  by  Rott  and  Lewellen  (45).  In  this 
work,  the  theories  for  both  rotating  and  stationary  disks  are  discussed. 

In  addition  to  the  review  of  past  research,  Rott  and  Lewellen  propose  an 
analytic  method  for  predicting  the  thickness  of  the  boundary  layer  and  the 
role  of  mass  ejection  for  the  case  of  a  turbulent  boundary  layer  and  a 
laminar  flow  field. 

The  theory  is  developed  by  assuming  universal,  turbulent,  boundary 
layer  velocity  profiles.  The  simplicity  of  the  method  makes  it  useful 
for  rapid  estimates  of  boundary  layer  effects  on  the  end-walls  of  a  rotating 
container.  As  in  other  previous  works  in  the  field,  the  velocity  of  the 
main  flow  field  Is  assuned  to  vary  in  the  following  manner. 

V  =  c  rn 


-  1  <  n  <  +  1 


The  results  of  this  theory  are  presented  In  Fig.  12.  It  may  be  noted 
that  when  a  free  vortex  exists,  the  radial  inflow  can  be  quite  large,  whereas 
when  n  f  -  1,  a  maxlmun  exists.  These  results  are  In  agreement  with  Mack 
(38,  39),  King  and  Lewellen  (40,  41,  42).  However,  these  results  differ 
from  the  results  of  Cooke  (33),  This  is  attributed  to  the  use  of  a  two 
parameter  boundary  layer  by  Cooke. 
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5.  BACKFLOWS  AND  VORTEX  BREAKDOWN 
5.1  Introduction 

The  phenomenon  of  backflow  or  reversed  axial  flow  in  a  vortex  has 
been  a  topic  of  interest  for  some  years.  A  classical  example  of  backflow 
is  found  in  the  Ranque-Hi lsch  vortex  tube.  In  this  case,  gases  are  In¬ 
jected  tangentially  into  an  open-ended  tube.  A  hot  gas,  rotating  in  one 
direction,  is  found  at  one  end,  while  a  cold  gas,  rotating  in  the  opposite 
sense.  Is  found  at  the  other  end.  This  effect  has  been  noted  in  a  number 
of  experiments  and  Is  analytically  predicted  with  the  models  of  Sullivan 
(22)  and  Donaldson  and  Sullivan  (16,  17,  18). 

Aside  from  straight  tube  flows,  backflows  have  been  noted  in  diverging 
and  converging  sections  for  both  initially  solid  and  free  vortices.  As 
yet,  the  analytical  description  of  the  process  is  incomplete  and  considerable 
controversy  exists  as  to  the  role  of  viscosity  in  setting  up  and  maintaining 
the  reverse  flow.  In  many  cases,  such  as  the  vortex  tube,  it  seems  that  the 
reverse  flow  Is  predominantly  a  viscous  effect.  There  are  indications,  how¬ 
ever,  that  the  concept  of  vortex  breakdown  In  which  a  small  stagnation 
region  of  locally  reversed  flow  near  the  axis  exists,  is  primarily  due  to 
dynamic  instability.  This  action  is  triggered  by  some  small  change  In  the 
boundary  conditions  so  that  the  flow  must  readjust. 

In  this  section,  the  concepts  of  vortex  breakdown  and  backflows  will 
be  discussed.  There  are  numerous  references  to  this  phenomena.  However, 
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only  a  few  of  these  will  be  presented.  Lord  Rayleigh  (46,  47,  48)  and  Lord 
Kelvin  (49)  wrote  about  the  conditions  of  stability  of  fluids  in  a  vortex 
flow  in  the  later  parts  of  the  19th  century.  The  primary  conclusion  of 
these  studies  is,  that  to  posses  dynamic  stability  In  the  presence  of  In¬ 
finitesimal  disturbances,  the  flow  must  have  a  circulation  which  constantly 
Increases  with  radius.  In  addition,  the  circulation  should  posses  the 
same  sign  along  a  radius.  When  this  Is  not  the  case,  instabilities  will 
be  developed,  as  in  the  well  known  case  of  rotating  coencentrlc  cylinders. 
This  concept  has  recently  been  expanded  by  Ludwelg  (50)  who  considers  the 
stability  of  an  Invlscld,  spiralling  flow  In  a  narrow  cylindrical  amnulus. 
This  study  Is  a  direct  extension  of  Rayleigh's  work  which  allows  radial 
gradients  to  exist  In  the  axial  velocity.  By  employing  small  disturbance 
theory,  the  stability  boundary  Is  deduced.  It  Is  shown  that  very  small 
radial  gradients  can  cause  a  flow  which  satisfies  Rayleigh's  criterion  to 
go  unstable.  The  work  of  Ludwleg  has  been  extended  by  Kiessling  (51)  to 
include  viscous  effects  for  a  flow  between  rotating  cylinders  of  small 
separation. 


5.2  Theory  of  Burgers  and  Weske 

In  two  articles  prepared  at  the  University  of  Maryland,  Weske  (52) 
and  Burgers  (53)  investigated  the  existence  of  backflow  as  derived  from 
"stretching"  a  rotational  vortex  core  in  converging  or  diverging  passages 
as  shown  In  Fig.  13.  The  analyses  presented  are  for  a  solid-body  vortex 
core  surrounded  by  an  irrotatlonal  vortex  which  flows  from  a  section  of 
uniform  radius  to  another  section  of  a  different  radius.  The  equations 
for  the  flow  are  derived  from  an  Incompressible,  Invlscld,  vortex  flow  In 
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a  cylindrical  coordinate  system . 

Weske  first  Investigated  the  problem  in  general  terms.  However, 
Burgers  expanded  the  analysis  and  presented  particular  solutions.  Burgers 
found  that  It  was  possible  to  achieve  exact  solutions  for  the  case  of  a 
core  with  solid-body  rotation  and  initially  uniform  axial  velocity  pro¬ 
files  surrounded  by  an  Irrotational  vortex.  The  solution  was  derived  from 
classical  Incompressible  fluid  mechanics  by  considering  the  following  set 
of  equations. 

LAMB'S  EQUATION 

V(P/p  +  1/2  q2)  =  q  x  u  (5.1) 

where, 

to  =  v  x  q  The  vorticity  vector  (5.2) 


CONTINUITY  EQUATION 


w 


1 

2  nr  dr  ’ 


_  djjj_ 

2ur  dr 


INITIAL  CONDITIONS 


w  =  W  ,  v  =  r  to. 


(5.3) 


(5.4) 


By  conserving  angular  momentum  along  streamlines  which  flow  from  one 
section  of  zero  radial  velocity  to  another,  the  following  equation  was 
derived. 


d2w. 


1 


dWr 


UL 


dr, 


2  +  2  .  +  W  ^w2  W) 


dr2 


(5.5) 


This  equation  has  a  Bessel  function  solution. 
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w2  =  W(1  -  a  i0  (  “g  %  )  (5.6) 

where , 

W  =  the  axial  velocity  at  the  Inlet 
o)0  -  the  angular  velocity  at  the  Inlet 
r2  *  the  radial  coordinate  at  the  second  straight  section 


The  constant,  a,  Is  an  undetermined  coefficient  which  can  be  evaluated 
from  conservation  of  mass  considerations. 


2ir 


5  Wr9  2ai  r 
w0r9dr5  =  ir(r2  +  a  J-|  (  jjj  )  ) 


2  2  2 
aW 


1  + 


wR*2  U1 


2uiR*p  p*  P 

'-it1)  -  <k> 


(5.7) 

(5.8) 


where, 


R*  =  the  Inlet  radius  of  the  rotational  core 
R*2  =  final  radius  of  the  core 


For  the  case  where  the  rotational  core  extends  throughout  the  flow  field, 
the  results  are  particularly  simple. 

2 


1-n 


not 


7"  TT^J 


w  =  w(i  +  (_ _ na  )  j 

2  wu  2  lF17TrtaF  Jo  [  W  } 
n  I 


(5.9) 

(5.10) 


n  =  R2/R 
a  -  2(i)qR/W 


where, 
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Thus,  It  can  be  seen  that  negative  axial  velocities  can  occur  under 
a  number  of  conditions. 


1. 

n2 

>  1 

2. 

Jo 

<  o  , 

J]  >  0 

3. 

J1 

<  0  , 

J1  >  0 

For  the  case  of  converging  flows,  n  is  less  than  one.  However,  JQ  and 
may  be  negative  depending  upon  their  arguments.  For  diverging  flows,  n  >  1. 
Backflow  can  be  expected  provided  the  Bessel  functions  remain  positive. 

In  addition.  It  Is  pertinent  to  add  that  no  negative  velocities  actually 
occur  unless  the  above  criteria  cause  negative  terms  which  are  larger  than 
the  initial  axial  velocity,  W.  It  Is  also  Interesting  to  note  that  multiple 
regions  of  back  flow  can  be  obtained  by  proper  choice  of  the  arguments  of 
the  Bessel  functions.  However,  these  become  successively  smaller  and 
finally  become  smaller  than  the  Initial  axial  velocity  term. 

The  more  general  case  of  the  rotational  core  surrounded  by  the 
Irrotational  vortex  Is  also  treated  In  some  detail  by  looking  at  limiting 
cases.  One  of  the  more  Interesting  cases  occurs  for  continuous  flow  with 
no  radius  change.  In  this  case,  the  Bessel  function  solution  is  appropriate 
for  radial  equilibrium,  however,  not  necessarily  stable  equilibrium.  The 
question  of  Instability  Is  studied  by  perturbing  the  Euler  equations.  It 
Is  found  that  the  perturbed  equations  allow  undamped  oscillations  to  occur. 
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5,3  Theories  of  Squire,  Brooke,  Benjamin  and  Lambourne 


In  1960,  Squire  (54)  proposed  a  mechanism  for  the  vortex  breakdown 
phenomenon.  It  was  postulated  that  when  the  flow  could  sustain  infinitesimal 
standing  waves,  a  characteristic  bubble  would  form  along  the  axis,  the  idea 
being  that  waves  generated  downstream,  or  at  the  exit,  would  collect  at  the 
critical  point.  However,  it  was  later  pointed  out  by  Benjamin  (55)  that 
the  group  velocity  of  such  waves  would  allow  only  upstream  disturbances  to 
be  felt  at  the  critical  point. 

In  1962,  Benjamin  (55)  published  his  theory  of  the  vortex  breakdown 
phenomena  which  was  closely  related  to  the  experimental  work  of  Harvey  (56). 
Harvey  had  described  the  formation  of  stagnation  regions  of  locally  reversed 
flow  on  the  center  line  of  a  vortex  tube.  The  theory  of  Benjamin  closely 
parallels  the  concept  of  the  hydraulic  jump  problem  where  a  flow  in  a  channel 
can  exist  in  two  conjugate  pairs,  namely  supercritical  and  subcritical.  The 
definition  of  the  critical  condition  is  that  condition  when  the  velocity  of 
the  long  waves  of  water  is  the  same  as  the  axial  velocity  of  the  flow  (Froude 
number  =  1.0).  Benjamin  hypothesized  that  an  analogous  condition  exists  in 
swirling  flow  such  that  there  can  be  a  sudden  transition  from  the  super¬ 
critical  to  the  subcritical  condition. 

The  theory  is  postulated  employing  the  calculus  of  variations  to 
establish  the  concept  of  conjugate  flows  based  upon  a  deficiency  of  the 
"force  flow",  S» 


S  = 


fr 

< 

o 


+  p)dA 
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Associated  with  the  conjugate  pair  of  flows  Is  a  gain  of  force  flow 
when  the  transition  from  supercritical  to  subcrltical  occurs.  Benjamin 
further  Investigated  the  conditions  for  the  onset  of  transition  which 
occurs  when  the  ratio  of  axial  to  tangential  velocity  is  of  the  order 
of  1 . 

A  subsequent  publication  by  Benjamin  (57)  was  presented  in  1964.  In 
this  paper,  the  essential  ideas  of  conjugate  flows  in  an  incompressible 
swirling  flow  are  reaffirmed.  In  1965,  Lambourne  (58)  published  an  analysis 
of  vortex  breakdown  phenomena.  The  theory  of  this  work  closely  parallels 
those  of  Weske  (52),  Burgers  (53),  and  Benjamin  (55).  In  this  paper  con¬ 
ditions  for  the  critical  conditions  are  derived  in  a  somewhat  different 
fashion.  However,  the  result  is  the  same  as  that  found  by  Squire  (54)  and 
Benjamin  (55).  The  analysis  proceeds  upon  the  basis  of  the  Bessel  function 
Solution  of  Burgers  (53),  In  this  study,  various  changes  were  made  in  the 
external  flow  and  regions  of  no  apparent  solution  were  found,  thus  suggesting 
a  finite  transition  or  vortex  breakdown. 

5.4  Theory  of  Gore  and  Ranz 

In  1 964,  Gore  and  Ranz  (59)  presented  a  paper  illustrating  backflow 
occurring  in  divergent  sections  of  a  flow  generated  by  a  rotating  porous 
plate.  A  theory  was  developed  employing  successive  approximations  to  a 
set  of  simplified  Navler-Stokes  equations.  Once  again  the  condition  of 
back-flow  occurred  when  the  ratio  of  the  axial  velocity  to  the  tangential 
velocity  at  the  wall  was  of  the  order  of  one.  Regions  of  back-flow  were 
predicted,  which  were  preceeded  by  a  stagnation  point.  This  phenomenon 
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has  direct  application  to  the  flame  pattern  of  a  swirling  flow  combustor, 

5,5  Conclusions 

In  conclusion,  the  existence  of  backflow  due  to  vortex  breakdown  has 
been  analytically  and  experimentally  verified.  The  theory,  however,  Is 
still  qualitative  since  the  phenomenon  is  related  to  the  stability  of  the 
flow.  The  characteristic  axisymmetric  bubble  of  reverse  flow  occurs  most 
frequently  in  the  diverging  section  of  a  supercritical  flow.  The  analytical 
and  experimental  work  of  Burgers  (53)  and  Weske  (52)  indicate  that  the 
phenomenon  may  also  occur  in  a  converging  section. 
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6.  BURNING  RATE  AUGMENTATION  IN  AN  ACCELERATION  ENVIRONMENT 

6. 1  Introduction 

The  rotation  imparted  to  a  solid  propellant  rocket  motor  not  only 
changes  the  flow  of  the  gases,  but  also  affects  the  burning  rate  of  the 
propel lanto  The  variations  in  burning  rate  have  been  studied  most  ex¬ 
tensively  by  experimental  methods  and  are  reported  upon  in  Ref.  69. 

However,  recently  a  few  analytical  investigations  in  this  field  have 
been  attempted  with  varying  success. 

6.2  Theory  of  Glick 

Glick  (60)  advanced  a  theory  to  account  for  the  effect  of  acceleration 
forces  on  the  burning  mechanism  of  solid  propellants.  This  was  accomplished 
by  extending  the  grannular  diffusion  model  of  Suranerfield  (61)  for  the 
burning  of  composite  solid  propellants  at  high  pressures.  Since  Sunmerfl eld's 
analysis  is  qualitative  and  empirical  in  nature,  the  extension  to  account 
for  acceleration  effects  should  also  be  viewed  in  the  same  manner. 

Two  basic  effects  are  considered  with  respect  to  acceleration  fields. 

1.  Pressure  gradients  across  the  gas  phase  reaction  zone. 

2.  Differential  forces  acting  upon  particles  of  different  density. 

The  first  possible  effect  was  eliminated  by  an  order  of  magnitude 

analysis.  It  was  shown  for  an  acceleration  field  of  50,000  g's  and  a  mean 
chamber  pressure  of  600  psia,  that  the  pressure  gradient  across  the  reaction 
zone  was  less  than  .01%.  Therefore,  only  the  effect  of  the  acceleration 
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field  acting  upon  the  density  inhomogenlties  was  incorporated  into  the 
analysis.  There  can  be  two  main  sources  of  density  inhomogenlties, 

1,  The  heterogeneous  nature  of  the  reaction  zone, 

2.  The  mean  temperature  gradient  through  the  reaction  zone. 

The  effect  of  these  density  inhomogenlties  due  to  the  acceleration 

forces  acting  normal  and  parallel  to  the  propellant  surface  were  investigated. 
It  was  concluded*  by  considering  the  theory  of  free  convection  flow  for 
the  case  of  normal  accelerations  and  boundary  layer  flow  for  the  case  of 
parallel  accelerations*  that  the  mean  temperature  gradient  has  only  secondary 
effect. 

Thus,  the  effect  of  acceleration  forces  is  primarily  concerned  with  the 
density  Inhomogenlties  produced  by  the  pockets  of  fuel  vapor  imbedded  in  a 
homogeneous  oxidizer  vapor.  In  the  high  pressure  limit,  the  granular  diffusion 
model  is  primarily  concerned  with  thickness  of  the  gas  phase  reaction  zone. 


where 

6r“  »  thickness  of  the  gas  phase  reaction  zone 
Vfv  -  normal  velocity  of  the  fuel  vapor 
Tfv  *  lifetime  of  the  fuel  vapor  pocket. 


(6.1) 


Since  the  burning  rate  Is  assumed  to  be  Inversely  proportional  to  the  thick¬ 
ness  of  the  gas  phase  reaction  zone,  the  effects  of  acceleration  on  burning 
should  be  explainable  in  terms  of  accelerations  effects  upon  Vf  and  Tfv* 
since  any  relative  motion  between  the  fuel  pocket  and  the  oxidizer  vapor 
will  Increase  the  interdiffusion  rate.  In  addition,  the  flow  of  the  oxidizer 
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vapors  can  be  reduced  If  the  acceleration  Is  Into  the  propellant.  Thus, 
two  effects  can  be  postulated:  a  scalar  effect  causing  an  increase  In 
Interdiffusion  which  is  independent  of  the  direction  of  the  acceleration 
field,  and  a  vector  effect  which  can  increase  or  decrease  the  thickness 
of  the  gas  phase  reaction  zone. 

Based  upon  these  considerations,  Glick  proposed  a  model  allowing  for 
both  a  vector  contribution  and  a  scalar  contribution  to  the  burning  rate. 

The  final  burning  rate  equation  was  written  in  terms  of  the  Reynolds 
number,  the  Grashof  nimber,  the  Schmidt  nunber,  and  the  diameter  of  the 
fuel  vapor  pocket. 

It  was  concluded  that  some  of  the  characteristics  noted  in  experimental 
observations  could  be  predicted.  However,  the  analysis  failed  at  high 
accelerations,  since  the  model  predicts  that  burning  rate  increased  with¬ 
out  bound  with  increasing  acceleration  whereas  experimental  results  in¬ 
dicated  a  maximum  burning  rate  with  g  level. 

6.2.2  Conclusions  from  G11ck‘s  Theory 

1.  "The  effects  of  acceleration  induced  pressure  difference  across 
the  mean  flow  of  the  gas  phase  reaction  zone  on  the  burning  rate 
of  nonmetallized  composite  propellants  are  negligible." 

2.  "The  major  effect  of  accelerations  on  the  burning  rate  of  a 
nonmetallized  composite  solid  propellant  is  derived  from  the 
effect  of  acceleration  on  the  heterogeneous  nature  of  the  gas 
phase  reaction  zone." 
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3„  "The  granular  diffusion  flame  model  correctly  predicts  the 

trends  for  the  variation  of  the  burning  rate  of  a  nonmetal llzed 
composite  propellant  with  respect  to  the  direction  of  the 
acceleration  vector  relative  to  the  burning  surface." 

4.  "The  granular  diffusion  flame  model  correctly  lead  to  an  upper 
bound  on  burning  rate  changes  when  the  acceleration  vector  Is 
normal  to  and  into  the  burning  surface." 

5.  "The  ratio  of  the  burning  rate  of  a  nonmetallized  composite 
propellant  In  an  acceleration  field  to  its  static  value  is  not 
dependent  upon  either  catalyst  content  or  initial  temperature 
when  the  acceleration  vector  is  parallel  to  the  burning  surface." 

6.  "The  ratio  of  the  burning  rate  of  a  nometal  llzed  composite 
propellant  In  an  acceleration  field  to  Its  static  value  Is 
generally  reduced  by  the  following:" 

a.  Addition  of  burning  rate  catalysts. 

b.  Reduction  of  the  mean  diameter  of  the  oxidizer  particles, 

c.  Increasing  the  initial  temperature  of  the  propellant. 

7.  "The  burning  rate  of  a  nometal  11  zed  composite  propellant  in  an 
acceleration  field  when  the  thickness  of  the  gas  phase  reaction 
zone  is  much  greater  than  the  diameter  of  a  pocket  of  fuel  vapor 
is  given  by 

1/2  ? 

„  G  ,  '  cose  ,  G  .  cos  e  c. 

f-- C,  -=4 - + 1  c,2  — — 2~  +  fr  3,/2 

°  cd  1/2  Re°  C“  Re°  0 


(6.2) 


Sh  =  the  Sherwood  number  -  hQ  d^/Dg 


where, 
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3  2 

Gf  s  a  •  Ap  dfy  /ug  =  The  Groshoff  nimber 
Re  -  the  Reynolds  nunber  =  p  p  rQ  dfy/yg 
2  “  constants  to  be  determined  from  experiment 
o  =  zero  acceleration 
a  -  acceleration 

C0  =  the  drag  coefficient  for  the  particle 

e  =  the  angle  of  the  acceleration  vector  to  the  propellant  surface. 

6, 3  Theory  of  Crowe  et  al 

In  a  final  report  for  a  contract  related  to  the  study  of  spin  stabilized 
rocket  motors,  Crowe  et  al .  (62)  discussed  various  models  for  burning  rate 
augmentation  due  to  motor  rotation.  These  are: 

1,  Effect  of  a  centrlbugal  force  field  on  the  gas  flow  from  a  solid 
propellant  burning  surface, 

2,  Contribution  of  particle  combustion  zone  to  surface  heating. 

3,  Retention  of  metal  oxide  particles  on  the  burning  surface, 

6.3,1  The  Effect  of  a  Centrifugal  Force  Field  on 
the  Gas  Flow  from  a  Solid  Propellant 

When  the  pressure  gradient  was  small,  the  momentim  and  energy  equations 
for  the  gas  flow  at  the  propellant  surface  in  an  acceleration  field  are 
given  by, 

u  -  a  =  centrifugal  acceleration  term  (6.3) 


82 


where , 


p  u 


dx 


u  + 

a 


d 

dx 


kd  h 


)  +  p  Q 


hQ  -  stagnation  enthalpy 
h  -  static  enthalpy 
k  thermal  conductivity 
Cp  -  specific  heat  at  constant  pressure. 


(6  =  4) 


By  examining  the  following  ratios,  some  estimate  of  the  effect  of  acceleration 
can  be  obtained. 


O'.S 

7 


where 


6  =  flame  stand  off  distance 
-  flame  temperature 

Tg  -  surface  temperature 
ion  g "s 


when 


S  <  1  irsrs 

0{a 


yw=c|! 


10“ 


(6.5) 

(6.6) 


(6.7) 

(6.8) 
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Thus,  no  significant  effects  of  rotation  should  be  felt  until  the  acceleration 
approaches  1000  g's.  Therefore,  it  was  concluded  that  the  acceleration  field 
does  not  drastically  affect  the  velocity  and  temperature  fields  at  the  pro¬ 
pellant  surface . 


6.3.2  The  Contribution  of  Particle  Combustion  Zone 
to  Surface  Heating 

While  the  acceleration  field  does  not  greatly  Influence  the  velocity 
and  thermal  fields  of  the  gases,  the  solid  particles  In  the  flow  may  be 
significantly  affected  by  acceleration.  When  the  acceleration  was  directed 
into  the  burning  surface,  the  metal  particles  were  retarded,  thus  enhaunclng 
heat  transfer  back  to  the  propellant  surface. 

The  distance  from  the  propellant  surface  to  any  particle  in  the  flow, 
assuming  Stokes  flow.  Is  given  by, 

S  =  (t  -  x  -  te"t/T)(u  -  at)  (6.9) 

where, 


t  -  time 


p  =  gas  viscosity 


It  may  be  noted  that  when, 


_  u 

T  -  — 

a 


(6.10) 


the  particles  will  remain  at  the  surface.  Thus  a  critical  particle  size 
can  be  deduced  for  retaining  the  particle  at  the  propellant  surface. 
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where 


pry  (1-w) 
r  -  [  §■  - ] 

P  2  Pp  °  Pg 


(6.11) 


w  =  aluminum  loading 


r  =  burning  rate 


ps  *  density  of  the  propellant 


In  order  to  determine  the  heat  transfer  at  the  surface,  the  combustion 
zone  was  broken  Into  four  regions. 


1.  Region  1  -  No  heat  addition  to  the  gas;  Q  =  0 

m  q 

2.  Region  2  -  Combustion  of  the  gas; 


3.  Region  3  -  Particle  Ignition; 

4.  Region  4  -  Particle  combustion; 


Q  = 


M 


q  =  o 
q  =  — r2- 


where, 


*Rg 


heat  of  reaction  for  the  gas 


qRp  =  heat  of  reaction  for  the  particles 

An  energy  equation  Is  then  postulated  assuming  the  kinetic  energy  to 
be  small  compared  to  the  thermal  energy. 


dfT  +  1 

V  k 


m  Cp  dT  _  1  dT 
k  dy  “  L  dy 


(6.12) 


This  equation  has  the  following  general  solution, 


T  =  C2J/l 


+.  L  S . J-  4.  r 

k  L1 


(6.13) 


where  and  C2  are  constants  of  integration. 
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The  above  equation  can  be  used  for  regions  2  and  4.  Combining  these 
solutions  and  employing  the  proper  boundary  conditions  for  each  region, 
the  temperature  at  any  point  can  be  determined  as  a  function  of  position,, 
It  was  concluded  that,  unless  the  solid  particles  are  in  the  order  of  one 
micron  or  less  above  the  surface,  no  appreciable  additional  heat  transfer 
will  occur  as  a  result  of  particle  combustion  zone  displacement. 


6.3.3  The  Retention  of  Metal  Particles  on  the 

Burning  Surface 


Based  upon  the  previous  analyses,  the  effects  of  the  acceleration 
field  must  be  primarily  in  the  retention  of  solid  particles  on  the  surface, 
since  only  negligible  effects  are  found  for  the  displacement  of  the  gas 
and  particle  combustion  zones. 

The  expression  for  a  critical  particle  size  as  derived  earlier  can 
be  written, 


_  .3  CD  Re  11  us,l/2 

rpc  '  (To  — ~ — > 

P 

For  Stokes  flow,  where  us  =  the  gas  velocity  at  the  surface, 


(6.14) 


pc 


2  Pp  a 


(6.15) 


Particles  larger  than  this  size  will  be  burned  at  the  surface  until 
they  reach  a  diameter  less  than  the  critical  diameter  whereupon  they  will 
move  into  the  gas  stream.  The  Increase  in  burning  rate  was  a  function  of 
the  additional  heat  released  to  the  propellant  surface  as  the  particle 


burns. 
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The  actual  burning  rate  law  was  postulated  as  a  ratio  of  the  burning 
rate  under  rotation  to  that  under  no  rotation. 


01 


1 


where. 


°  w_,  l  r  ^ 

«=0  1  -  -E—  H(-^  .  ». 

nr  pm  pm 


(6.16) 


H  =  a  function  to  be  determined 
hf  =  the  heat  of  vaporization  of  the  particles 
wm  s  the  propellant  metal  loading  weight  fraction 
L  =  the  heat  transfered  per  unit  mass  of  particle  burned 
rpm  =  mass  med'an  Particle  diameter 
*pm  s  Particle  burning  time  of  a  mass  median  particle 
o  ~  the  standard  deviation  of  the  particle  size  distribution. 


This  equation  is  used  as  a  tool  for  extrapolating  existing  data, 
since  the  constants  required  and  the  function  H  are  not  generally  known. 


6.3.4  Conclusions  from  Crowe's  Theory 

In  conclusion,  the  main  point  of  this  analysis  was  that  the  major 
effect  of  the  acceleration  field  was  to  move  particles  larger  than  the 
critical  size  to  the  propellant  surface.  The  particles  at  the  surface 
burn  and  release  their  energy  at  the  propellant  surface.  The  residence 
time  of  the  particle  at  the  surface  was  a  function  of  the  heat  release 
rate,  the  particle  size  and  the  acceleration  level.  Further,  It  was 
predicted  that  when  the  mass  median  particle  diameter  was  lower  than  the 
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critical  diameter,  the  effects  of  acceleration  will  be  small.  In  addition, 
if  there  are  no  metal  oxide  particles,  there  will  be  no  effect  regardless 
of  acceleration  level. 


6,4  Conclusions 

The  two  theories  reviewed  here  are  attempts  to  get  at  the  most  Important 
factors  Influencing  the  burning  rate  augmentation  in  an  acceleration  en¬ 
vironment.  Gllck's  theory  Is  primarily  for  nonalumized  propellants  and 
consists  of  a  scalor  and  a  vector  effect.  The  scalar  contribution  was  due 
to  the  Increased  relative  motion  between  the  fuel  rockets  and  oxidizer  vapor. 
The  vector  effect  was  associated  with  the  compression  or  expansion  of  the 
gas  phase  reaction  zone.  The  theory  of  Crowe  et  al.  was  based  entirely 
upon  metalized  propellants  and  the  effect  of  retaining  larger  metal  particles 
upon  the  propellant  surface.  A  further  description  of  the  validity  of  these 
theories  in  comparison  to  experimental  results  Is  given  in  a  theses  by 
Anderson  (63)  and  Northam  (64)  which  are  reviewed  with  many  other  in  reference 
(69). 
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7 .  SURVEY  ARTICLES 
7.1  Swithenbank  and  Sotter 

In  1964  Swithenbank  and  Sotter  (65,  66)  presented  an  analysis  of  the 
effect  of  vortlcles  generated  by  combustion  pressure  oscillations  on  the 
performance  of  solid  propellant  rocket  motors*  The  authors  Investigated 
the  theories  of  vortex  flow  to  determine  their  applicability  in  the  chamber 
of  a  rocket  motor,  and  they  also  investigated  Mager's  (3)  analysis  for 
rotating  nozzle  flows* 

Applying  incompressible  viscous  flow  theory  in  the  chamber,  an  estimate 
was  made  of  the  viscous  effects  upon  the  vortex  distribution.  For  a  typical 
six  Inch  diameter  rocket  motor,  the  turbulent  Reynolds  nimbers  for  the  radial 
and  tangential  velocities  were  estimated.  Based  upon  the  turbulent  velocity 
correlation  of  Keyes  (67),  these  Reynold's  nunber  were  estimated  to  be 

Re>t=2xl05  Re»r  =  8,0  (7J) 

Thus,  It  was  concluded  that  the  vortex  motion  will  not  appreciably 
be  affected  by  viscous  forces  except  possibly  near  the  centerline.  Employing 
the  incompressible  flow  theory  of  Einstein  and  Li  (15),  it  was  found  that 
the  flow  will  be  predominately  free  vortex  in  character  except  near  the 
centerline  where  a  solid  body  vortex  should  exist. 

For  the  nozzle  flow*  viscous  Incompressible  flow  theory  was  discarded 
In  favor  of  the  invlscid  compressible  flow  analysis  of  Mager  (3).  That 
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theory  predicted  a  decrease  In  the  mass  flow  and  thrust  producing  capabilities 
of  a  nozzle  due  to  the  following  factors: 

1.  A  vacuum  core,  region  on  the  centerline,  and 

2.  A  reduction  In  the  mass  flux  at  the  throat. 

For  the  measurements  of  the  vortlcles  generated  In  their  experimental 
studies,  Swlthenbank  and  Sotter  found  the  mass  flow  to  be  reduced  70%  when 
Mager's  theory  Is  applied.  In  addition  to  predicting  the  reduction  of 
effective  throat  area,  Mager's  theory  was  used  to  predict  observed  Irregular 
burning  caused  by  vortex  flow.  A  decrease  in  effective  throat  area  can  be 
related  to  the  chamber  pressure  in  the  following  manner. 


where, 

Ay  =  area  of  the  burning  surface 
r  =  Pn  =  the  burning  rate 
Pp  -  density  of  the  propellant 
An**  *  the  effective  throat  area 
Pc  -  the  chamber  pressure 
v  =  characteristic  velocity  of  the  gases. 


(7.2) 

to 


The  following  conclusions  regarding  the  affects  of  rotation  were  made: 

1.  The  decrease  in  effective  throat  area  increased  the  chamber 
pressure. 

2.  The  burning  rate  was  Increased  due  to  increased  chamber  pressure. 

3.  Erosive  effects  of  the  vortex  further  increased  the  burning  rate. 
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7,2  Crowe 

A  final  report  by  Crowe  (62)  released  In  1966  presents  an  excellent 
study  of  the  many  theories  discussed  in  this  literature  survey.  In  addition, 
an  Investigation  was  made  o *  the  dynamics  of  particle  flow  In  vortic1es„ 

The  experimental  portion  of  this  research  program  Is  reviewed  in  Refs,  69, 
Therefore,  only  the  analytical  investigation  will  be  discussed  here. 

The  analysis  of  the  flow  in  a  rocket  motor  was  divided  into  the  following 
sections: 

TL  Solution  for  the  viscous  vortex  motion  in  the  rocket  chamber, 

2,  Distribution  of  the  metal  oxide  particles  in  the  rocket  chamber, 

3,  Effect  of  rotation  upon  nozzle  performance. 


7,2,1  Solution  for  the  Viscous  Vortex  Motion 
In  the  Rocket  Chamber 


The  analysis  of  the  viscous  vortex  in  the  chamber  was  prefeced  by  a 
solution  for  the  fnviseld,  compressible  flow  case.  The  Invlscld  compressible 
flow  case  was  used  as  a  eonspaHson  to  estimate  the  effects  of  compressibility. 
The  Euler  equations  are  presented  ana  solved  in  general  terms  under  the  con¬ 
ditions  of  conservation  of  angular  momentum »  assuming  that  all  gas  particles 
possess  the  same  angular  momentum  since  they  originate  at  the  same  radial 
position. 


y 


v  - 


V  T 

Vi 


(7,3) 


w 


at  the  propellants  surface 
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r  =  radius  of  the  propellant  surface 

W 

Since  the  motion  of  the  vortex  Is  irrotatlonal  and  the  flow  is  invlscld 
and  adlabadlc,  the  flow  properties  are  governed  by  an  isentropic  expansion 
from  the  Initial  conditions- 

The  maxlmwi  velocity  occurs  when  the  static  temperature  becomes  zero- 


w  ■  <7-5> 

y-1 

When  the  flow  Is  viscous  in  nature,  the  shear  stress  will  alter  the 
free  vortex  solution-  The  shear  stresses  in  the  tangential  direction  are 
given  by 


/3V  V\  =  0 

T  ~  y  (aF  “  F  ~  "  2  u 


v  r 
w  w 


(7.6) 


Thus  as  the  radius  decreases,  the  shear  stresses  increase. 


t  a  (1/r)2  .  (7.7) 


The  viscous  torque  per  unit  length  parallel  to  the  flow  becomes 

M  -  +  (zTtr)x  °  r  =  -  4 ivy  v  r  .  (7.8) 

W  W 

Thus,  the  potential  vortex  gives  a  constant  torque  distribution.  Because 
there  existed  a  free  boundary  near  the  centerline,  the  torque  at  that  point 
must  vanish.  Therefore,  a  solid  body  vortex  should  have  existed  near  the 
centerline. 

The  approach  to  the  viscous  flow  solution  was  that  given  by  Burgers  (11). 
In  that  case,  the  axial  velocity  was  of  the  form 


W  =  -  2  uw  z/rw 


(7.9) 
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From  the  continuity  equation 


u  =  <r/rw>  “w  • 


(7.10) 


This  procedure  represents  the  stagnation  vortex  approach.  Following  the 
method  of  Burgers  (11)  and  Rott  (13,  14),  which  is  a  special  case  of  the 
Donaldson  and  Sullivan  solutions  (16,  17,  18,  22),  the  tangential  velocity 
was  given  by 


1  - 

1  -  e-Re/4 


(7.11) 


where 


n  =  r/r, 


w 


Re 


uw 


2  r. 


w 


Pressure  variations  for  the  Incompressible,  viscous,  vortex  flow  field 
in  the  chamber  are  obtained  by  employing  Rott's  analysis  (13,  14),  which  is 
an  extension  to  Burger's  solution  (13,  14). 


P  =  P  + 
c 


where 


A 


[i 


4(f)2] 

rw 


rl  ..2 


dn 


(7.12) 


Pc  =  pressure  at  the  wall  and  head  end  of  the  grain. 


The  term  on  the  right  hand  side  Involving  the  integral  may  be  evaluated 
employing  the  expression  derived  earlier  for  the  tangential  velocity. 
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P  ,  pc  ♦  p  -  .  t,  -  u2  -  4  <f-)2]  (f)  f(f^  .  n) 

w 

,2 

Re/4  2 


f(y) 


-Re/4 


Lhl 


w  '  (1-e‘Re/4)‘ 

■Re/4  -e~r  ,,  -yS  1  4 

C2-^y]  +  V-fe 


(7.13) 


+  [El  (y2)  -  El  '#)]  '-  [El  ( 2y2)  -  El  (Re/2)] 


where 


Ei(x)  *  I  “  dt, 

*  n  ‘ 


*^e 


(7.14) 


In  conclusion,  the  tangential  velocity  and  the  pressure  functions  were 

6 

plotted  for  various  values  of  Reynolds  numbers  from  1  to  10  . 

The  higher  the  Reynolds  number,  the  more  the  flow  approaches  a  free 
vortex.  For  Reynolds  numbers  of  the  order  of  one,  the  entire  flow  field 
becomes  forced  vortex  In  character.  Typical  values  of  Reynolds  nunber 
were  postulated  to  be  from  104  to  10*\  thus  predicting  that  the  tangential 
velocity  was  primarily  free  vortex  In  character,  while  the  pressure  gradient 
no  longer  Increases  with  Reynolds  number,  as  Illustrated  In  Figures  14  and 
15. 

It  should  be  noted  that  the  Reynolds  nunber  was  based  upon  the  radial 
velocity  at  the  propellant  surface.  The  radial  velocity  was  usually  In 
the  order  of  a  few  feet  per  second.  In  addition,  the  effects  of  turbulence 
are  to  increase  the  effective  viscosity  so  that  values  of  the  radial  Reynolds 
number  are  in  the  order  of  10  (65,  67).  Thus  in  an  actual  flow,  the  solid 
body  core  may  be  a  considerable  portion  of  the  total  flow. 
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FIGURE  14  THE  TANGENTIAL  VELOCITY  DISTRIBUTIONS  AS  A  FUNCTION  OF 
REYNOLD'S  NUMBER  USING  BURGERS  SOLUTION  APPLIED  TO 
INTERNAL  BURNING  MOTORS  (02). 

(Taken  from  reference  (62)  ,  Crowe  et  al.) 


7.2.2  Distribution  of  the  Metal  Oxide  Particle 
fn  the  Rocket  Chamber 


The  treatment  of  the  metal  oxide  particles  in  the  chamber  was  based 
upon  the  following  assumptions: 

1.  The  gas  properties  are  Independent  of  the  particle  pro¬ 
perties  and  are  those  determined  from  a  viscous,  Incompresslbl 
flow  analysis. 

2.  The  particles  are  In  angular  momentum  equilibrium  with  the 
gasses. 

3.  The  particles  do  not  penetrate  the  viscous  portion  of  the 
flow. 

4.  Stokes  flow  Is  assumed  for  the  particles. 

The  following  equations  of  motion  for  the  particles  were  derived. 


du  u  1  /  \ 

dt  "  nr  =  T  u6Wn  “ 

W 


dv  ,  uv 

”  at  +  — 

ID 


-  (cd  r 
T  W 


v  ) 

n 


(7.15) 

(7.16) 


~  [2  ugw  +  id]  . 
w 


(7.17) 


Solutions  for  these  equations  were  obtained  analytically,  and  the  distri¬ 
bution  of  the  metal  oxide  particles  was  found.  The  results  were  plotted 
for  typical  particle  diameters. 

A  critical  radius,  n_,  was  found  below  which  the  solid  particles  can- 

V 

not  penetrate  due  to  the  centrifugal  forces.  The  critical  radius  was  pro¬ 
portional  to  the  particle  size  and  the  square  root  of  the  motor  spin  rate. 
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For  parttcles  of  the  order  of  one  micron,  the  maldistribution  of  particle 
mass  flux  was  small  for  a  spin  rate  of  1000  rpm  or  less,  but  could  become 
appreciable  for  speeds  greater  than  10,000  rpm. 


7.2.3  Effects  of  Rotation  on  Nozzle  Performance 


The  analysis  of  the  nozzle  flow  condition  was  based  primarily  upon 
Mager's  analysis  (3)  of  Irrotatlonal  vortlcles  In  rocket  nozzles.  There¬ 
fore,  the  rotational  core  predicted  for  viscous  Incompressible  flow  was 
neglected  by  assuming  that  It  has  a  small  area.  Since  the  particles  In 
the  flow  do  not  pass  Into  the  viscous  region,  the  particle  distribution 
can  be  found. 

Since  the  particle  distribution  was  determined  at  the  nozzle  Inlet, 
it  was  only  necessary  to  apply  the  radial  momentum  equation  for  the  particles 
to  determine  their  redistribution  in  the  nozzle.  The  change  In  radial  posi¬ 
tion  Is  due  to  change  in  the  centrifugal  forces  acting  upon  the  particles, 
under  the  assumption  of  angular  momentum  equilibrium.  Thus 


CD  Re 


P  V 

MP  P 


2  <ug 


(7.18) 


where,  the  first  term  on  the  right  hand  side  is  a  centrifugal  force  term 
and  the  last  term  Is  the  Stokes  flow  viscous  term  of  the  particle  velocity 
lag  in  the  radial  direction. 

In  conclusion  It  was  found  that  when, 

Rt  =  1.0  inches 
Pc  =  1000  psia 


dp  =  1  micron 
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the  effect  of  a  rotating  flow  in  a  nozzle  upon  the  particle  flow  is  not 
large  for  spin  rates  below  6500  rpm.  When  the  spin  rate  was  larger,  the 
effect  on  the  particles  Increases.  In  addition,  when  the  particle  diameter 
was  above  1  micron,  the  effect  was  accentuated  so  that  an  appreciable 
fraction  of  the  particles  are  located  at  or  near  the  wall. 
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8.0  SUMMARY 

The  field  of  contained  rotating  flows  is  rapidly  growing  and  quite 
diverse.  In  this  report  only  some  of  the  possible  areas  of  study  have 
been  examined.  The  primary  emphasis  has  been  placed  upon  the  problem 
areas  which  would  be  found  in  rotating  rocket  motors.  It  is  evident 
from  the  number  and  variety  of  the  articles  listed  that  a  complete  des¬ 
cription  of  the  flow  in  the  nozzle  and  chamber  of  a  rotating  rocket 
motor  is  not  yet  possible.  It  should  be  noted  that  the  analytical  con¬ 
veniences  of  infinite  cylinders  and  rotating  plates  do  not  touch  upon 
the  Interaction  of  the  nozzle  and  the  chamber  flow  fields. 

In  general  the  effect  of  rotation  is  to  reduce  the  mass  flux  at 
the  throat.  This  effect  is  accentuated  by  increasing  the  motor  spin 
rate  and  by  increasing  the  nozzle  contraction  ratio.  However,  a  signifi¬ 
cant  portion  of  this  decrease  may  be  negated  by  a  transfer  of  angular 
momentum  to  the  nozzle  walls  through  viscous  shear. 

It  appears  possible  that  recirculating  viscous  flow  patterns  can  be 
established  from  the  action  of  the  fore  and  aft  end  sections  of  the  motor. 
Also,  the  boundary  layer  can  transport  significant  portions  of  the  mass 
from  one  area  of  the  flow  field  to  another.  There  Is  also  a  case  for 
back  flow  patterns  due  to  dynamic  effects  in  the  nozzle.  These  patterns, 
as  predicted  analytically,  can  be  quite  complex  if  the  ratio  of  tangential 
to  axial  velocity  becomes  much  greater  than  unity.  Portions  of  these 
patterns  tend,  however,  to  be  unstable  to  distrubances. 


100 


The  effect  of  the  reduction  In  mass  flow  rate  at  the  throat  Is  felt 
by  the  chamber  as  an  effective  decrease  In  throat  area.  This  can  cause 
Increases  In  chamber  pressure  over  the  non-rotating  condition.  This  In¬ 
crease  In  chamber  pressure  can  be  quite  severe  If  the  sensitivity  of  the 
propellant  to  pressure  Is  large.  The  effect  of  backflow  and  reclrculatory 
patterns  Is  to  produce  variations  in  the  local  burning  rates  where  these 
hot  combustion  gases  tend  to  impinge  upon  the  propellant. 

In  the  chamber  the  effect  of  rotating  flow  Is  to  set  up  pressure 
gradients  which  may  alter  the  local  burning  rates.  In  addition,  the 
acceleration  forces  accompanying  the  pressure  gradients  can  cause  retention 
of  solid  phases  at  the  propellant  surfaces  or  otherwise  affect  the  com¬ 
bustion  zone  of  the  propellant. 
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NOMENCLATURE 


English  Symbols 
a  -  accoustlc  speed 

aQ  -  stagnation  accoustlc  speed 

C  -  V  »  r  “  angular  momentum 

cp  -  specific  heat  at  constant  pressure 

D  -  diameter 

H  -  enthalpy 

m  -  mass  flow  rate 

n  -  burning  rate  Index;  tangential  velocity  exponent  v  =  k/rn„ 

P  s  pressure 

Q  s  volume  flow  rate 

q  s  the  total  velocity 

R  s  gas  constant 


\ 

Re 

r 

S 

T 

u 

V 

w 

z 


-  radius  of  the  wall 

-  Reynolds  nun  be r 

s  radial  space  coordinate 

-  entropy 

-  temperature 

-  radial  velocity 

-  tangential  velocity 

-  axial  velocity 

-  the  axial  space  coordinate 


1 02 


Greek  Symbols 

r  -  angular  momentim 

y  s  specific  heat  ratio 

6, a  =  boundary  layer  thicknesses 

g  ~  eddy  viscosity 


c 

n 

e 

u 

V 

IT 


-  ifi  reduced  axial  coordinate 

=  (r/r0)  reduced  radial  coordinate 
=  the  tangential  coordinate;  flow  angle 

-  absolute  viscosity 

~  u/p »  kinematic  viscosity 

-  3.14159 


P 

a 


♦ 

u> 

u 


-  density 

=  flow  ratio  -  u  r  p/u 

-  shear  stress 

”  velocity  potential  function 

-  stream  function 

-  angular  velocity 

-  vortlclty  vector 


Subscripts 

1  *  initial  condition 

o  -  stagnation  condition;  total 

p  -  particle 

x  -  tangential;  total 

w  -  wall 

eff  "  effective 

*  “  reduced  conditions  throat  condition 
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